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Abstract 

In a stationary axisymmetric spacetime, the angular velocity of a station¬ 
ary observer that Fermi-Walker transports its acceleration vector is also the 
angular velocity that locally extremizes the magnitude of the acceleration of 
such an observer, and conversely if the spacetime is also symmetric under 
reversing both t and (p together. Thus a congruence of Nonrotating Accelera¬ 
tion Worldlines (NAW) is equivalent to a Stationary Congruence Accelerating 
Locally Extremely (SCALE). These congruences are dehned completely lo¬ 
cally, unlike the case of Zero Angular Momentum Observers (ZAMOs), which 
requires knowledge around a symmetry axis. The SCALE subcase of a Sta¬ 
tionary Congruence Accelerating Maximally (SCAM) is made up of stationary 
worldlines that may be considered to be locally most nearly at rest in a station¬ 
ary axisymmetric gravitational field. Formulas for the angular velocity and 
other properties of the SCALEs are given explicitly on a generalization of an 
equatorial plane, infinitesimally near a symmetry axis, and in a slowly rotat¬ 
ing gravitational field, including the far-held limit, where the SCAM is shown 
to be counterrotating relative to inhnity. These formulas are evaluated in par¬ 
ticular detail for the Kerr-Newman metric. Various other congruences are also 
dehned, such as a Stationary Congruence Rotating at Minimum (SCRAM), 
and Stationary Worldlines Accelerating Radially Maximally (SWARM), both 
of which coincide with a SCAM on an equatorial plane of rehection symmetry. 
Applications are also made to the gravitational helds of maximally rotating 
stars, the Sun, and the Solar System. 
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1 Introduction 


A stationary axisymmetric spacetime has various preferred congruences of stationary 
observers, such as those whose four-velocities are parallel to the Killing vector held 
that is timelike at radial inhnity (wordlines nonrotating relative to inhnity), or 
those whose four-velocities are perpendicular to the Killing vector held which has 
closed orbits and which vanishes on the symmetry axis (Zero Angular Momentum 
Observers, or ZAMOs) BB- Here a new preferred congruence is dehned (SCAM, 
a special case of SCALE = NAW) in terms of the purely local properties of the 
commuting Killing vector helds, without reference to what they do elsewhere (e.g., 
at radial inhnity or around the symmetry axis). 

Using the MTW sign conventions 0 — in particular, the metric sign convention 
(--I--I--I-) — and the same boldface symbols for vectors and for the corresponding 
one-forms that have components obtained by using the metric tensor to lower the 
vector components, consider a region of spacetime with two independent Killing 
vector helds, vector helds 


k = k^djdx^, 1 = rO/dx^, 

(1) 

that are independent (not obeying ak -|- 61 = 0 for any 
zero) and whose corresponding 1-form components 

constants a and b not both 

9a(3^ 5 Qaf}^ j 

(2) 

obey Killing’s equation. 


^a;/3 ^/9;Q!) ^a;l3 ^/3;a- 

(3) 


Assume that these two Killing vector helds k and 1 also have the following three 
additional properties (though only the hrst two properties are necessary for the hrst 
part of the theorem to be proved): 

(1) The 2-form 

A = k A 1 (4) 

corresponding to the Killing bivector is timelike, obeying 

= 2(k-k)(M) - 2(k-l)2 < 0. (5) 

Then in each sufficiently small neighborhood one can redehne, if necessary, k and 1 
to be two new independent linear combinations of the original Killing vectors such 
that 1 is spacelike and the orthogonal vector held (lT)k — (kT)l, which by Eq. (^) 
is necessarily timelike, is future pointing (by choosing the appropriate sign for k). 
Then the linear combination 

K = k + m (6) 
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is future-pointing timelike at each point in the neighborhood for some hnite range 
of the constant say f 2 _i < < fl+i- 

(2) The two Killing vector helds commnte, 

[k, 1| = 0, (7) 

or, in component form, 

This implies that one can choose two of the fonr coordinates, say 
snch that 

k = d/dt, 1 = djdip. 

(3) The 2-form A = k A 1 obeys 

* A A =t:dA = 0, (10) 

which, when A 7 ^ 0 holds, as is implied by Eq. (|^), is eqnivalent to the orthogonally 
transitive or circularity condition 

kAlAdk = kAlAdl = 0. (11) 

In component form it is eqnivalent to 

= 0 ( 12 ) 


( 8 ) 

= t and = Lp, 


(9) 


or 


pk^-^g\ l[oikf}l'y\S\ 0 - 


(13) 


This condition is eqnivalent to the condition that one may constrnct in a local 
neighborhood a family of two-snrfaces orthogonal to both Killing vector helds P, 

1^. One may dehne two coordinates, say = (r, 6) for a = 2, 3, on these two-snrfaces 
such that orbits of the Killing vectors k and 1, and hence of all stationary observers, 
each stay at hxed x°‘. Then d/dr and d/dd are both orthogonal to k = d/dt and to 
1 = d/dip^ so in this coordinate basis the metric tensor has no components mixing 
the hrst two (0 or 1 ) and the last two (2 or 3) indices, he., it is block diagonal. As 
a result, A may be written as 


A = —Ddt A d(p. 


(14) 


where 

- D = go^gii - goigiQ = gug^^ - gt^g^t = = (k-k)(M) - (kd)^ < 0 (15) 

is the determinant of the hrst two-dimensional block of the metric. 
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Another simple way to state this third condition is to say Q that the spacetime is 
invariant under the simultaneous reversal of both coordinates t and ip. This follows 
from the block diagonality of the metric, and it implies that each of the quantities 
in Eq. ([TT|) are zero, since they are odd under this transformation. 

The most important examples of spacetimes with two independent commuting 
Killing vectors obeying these two properties are asymptotically flat stationary ax- 
isymmetric spacetimes with the Ricci tensor obeying 


k^R^lakf^lp] 0 , 


(16) 


which implies that property (3) above holds [§, P, |T3, though Eq. (|T^ just by itself 
does not imply property (3). In such spacetimes one may uniquely choose the Killing 
vector helds such that k = d/dt is & unit timelike vector held at radial inhnity and 
1 = d/d(f is a spacelike vector that vanishes on the symmetry axes (e.g., at 0 = 0 or 
6 = Ti) and has closed orbits with period Aip = 27r. However, the results below apply 
more generally, assuming only that k and 1 have k A 1 timelike and obey Killing’s 
Eq. (^, the commutativity condition [k, 1] = 0, and the two-surface-orthogonality 
condition =t:(k A 1) A *d(k A 1) = 0. 

An observer whose four-velocity is 


u = 


(17) 


with K = k -|- HI with hxed H = d<p/dt (which shall be called the angular velocity, 
since that is what it for a stationary axisymmetric spacetime) may be dehned to be 
a stationary observer (SO). A stationary congruence of observers (SCO) is a space- 
hlling family of observers (one crossing each point of each local spatial hypersurface 
in the region of spacetime under consideration) with four-velocities 


(k + H(3;^)l) 
|k -I- H(a;“)l| 


(18) 


that have the angular velocity H depending only on the two coordinates that stay 
hxed along each worldline. 

In order to calculate the acceleration vector 


a = Vuu 


(19) 


with H constant along the worldline, it is convenient to consider the Killing vector 
held K = k -|- HI with this same H hxed as a constant everywhere (and not having 
the spatial variation with x°‘ that a stationary congruence of diherent worldlines at 
diherent x°‘ might have). If for this hxed (^^-independent) H one dehnes the scalar 


held 


$ = $(H) = ^ln(-K-K), 


( 20 ) 
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then 


u = e '*’K 


( 21 ) 


is, over the region of spacetime where K is (future-pointing) timelike, the four- 
velocity of a rigidly rotating stationary congruence of observers, differing from the 
four-velocities U (7 of the congruence given by Eq. (^), where hi is allowed to be 
a function of (i.e., different angular velocities for different stationary observers 
within the congruence, though I am always taking the angular velocity hi to be hxed 
for a given observer). 

Then the antisymmetry of the covariant derivative of the Killing vector field K 
implies that the covariant components of the acceleration vector are 


a„ = u^Ua-,p = e '*’i^a);/3 

= e-^^[~{K^K0),^ + K^KPK^K^,0/{-K-K,)] 

= ( 22 ) 


which are nonzero only for a = a = 2 or 3. That is. 


a = V<h (23) 

is perpendicular to both k and 1. This fact requires property (2) but not property 
(3) above. 


2 Stationary Congruence Accelerating 
Locally Extremely (SCALE) 


Now for each value of the pair of coordinates x“, we would like to hnd the value of 
n that extremizes the magnitude of the acceleration of the corresponding stationary 
observer. For this purpose, it is convenient to define the following scalar fields 
(functions of a:“): 


F = FiSl) 


A = -k-k = -goo = -gtt, 

(24) 

B = -k-l = -goi = -gt^, 

(25) 

C = IT 5^11 = 

(26) 

= e^'*’ = -K-K = 21 + 2VLB + 

(27) 
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( 28 ) 


G = G(Si) = ~F(Si) = -K-I = B + aC, 

2 c/iZ 

where the partial derivative with respect to is at hxed and so at hxed A, B, and 
C. Since we have chosen 1 to be spacelike, C < 0. If the source of the axisymmetric 
gravitational held is entirely rotating in the same direction, B will typically have 
the sign of this direction everywhere, which by the appropriate choice of the sign 
of the coordinate tf can be chosen to be positive. A will be positive outside any 
ergospheres but negative inside, if any exist. However, property (2) above implies 
that we are outside any Killing horizon, so 

D = B^-AC = —A^^Ao^p > 0. (29) 


Then 


so 


dsL 


V(f) 


and thus 


$ = 



(30) 


^ VF VH + 2HVF + 

2{A + 2nB + n‘^c) ’ 

{AVB - BVA) + H(HVC - CVA) + ^^{BVC - CVB) 
{A + 2QB + H2C')2 


(31) 

(32) 


da^ 

'm 


(9a 1 C 

2a-— = -VF-V(-) = F-^VF-{FVG-GVF) 
oil F F 

2V<I>- V(K-1/K-K) = (K-K)-V(K-K)-V(K-1/K-K) 

[VA+2QVB+Q‘^VG] ■ [{AVB-BVA)+Q{AVG-GVA)+Q\BVG-GVB)] 

{A + 2QB + Q'^Gf 


(33) 


using [K, 1] =0. 

Setting da^/dfl to zero at each point gives an extremum of the magnitude of the 
acceleration of a stationary observer there, with H thus obeying the equation 

[VA+2QVB+Q^VG]-[{AVB-BVA)+Q{AVG-GVA)+Q\BVG-GVB)]=0. (34) 

Expanded out in powers of the angular velocity H, this is a quartic equation for 
H(x“), with coefficients that are combinations of H = —goo, B = —goi, G = —gn, 
and dot products of their gradients. 

A congruence of stationary worldlines corresponding to one of the roots of Eq. 
(|34[) for f2(a;“) at each might be called a Stationary Congruence Accelerating 
Locally Extremely (SCALE), and if the local extremum of the acceleration (as a 
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function of f2) is a (local) maximum, the congruence might be called a Stationary 
Congruence Accelerating Maximally (SCAM). Since in the frame of the device (e.g., 
a rocket) accelerating the observer, the magnitude of the acceleration may be inter¬ 
preted as the apparent weight or heaviness of the observer (e.g., as in saying that an 
observer in free fall along a geodesic is “weightless”), one might say that an observer 
moving along a SCALE is “extremely heavy,” taking extreme to mean either a local 
maximum (for a SCAM) or a minimum (for the other SCALEs). 

(Generically there is no global maximum for the acceleration, since it can be 
made arbitrarily large by making arbitrarily near one of the two endpoints 
and fl+i of its allowed range. 


— 


B ± - AC 


(35) 


where F goes to zero and hence K = k -f fil becomes null, unless this endpoint 
corresponds to a null geodesic where VE also goes to zero, in which case the ac¬ 
celeration stays hnite. This last fact uses property (1) in the form of Eq. (^^ , 
D = — AC > 0, so that the two endpoints have a nonzero separation in and 

hence F, as a quadratic polynomial of given by Eq. (pTl), has only a simple zero 


at each end and cannot give inhnity when divided into VF if the latter also has a 
zero at the corresponding endpoint.) 

Although typically a? thus has no global maximum within the allowed range of 
where K is timelike, there are usually (at least in weak gravitational helds) two 
local minima for a? and one local maximum between these two minima, though it is 
also possible in strong gravitational helds to have only one local minimum and no 
local maximum. 

For a stationary axisymmetric spacetime in a region of weak gravity, one of the 
roots of the quartic Eq. (Q) in Q{x°') corresponds to an imaginary spatial four- 
velocity u = i(K-K)“^/^K (unphysical), and the other three roots correspond to 
real timelike four-velocities, with the two outer roots (say and f2+) giving local 
minima of the acceleration and the root in between (say flo) giving a local max¬ 
imum. In a static spacetime, the local maximum of the acceleration occurs for a 
static worldline, at flo = 0, accelerating against the pull of gravity to stay at a 
hxed position. In a Newtonian description in which 1 = d/d(p vanishes along the 
z-axis, the two local minima of the acceleration occur at the angular velocities 
at which the centrifugal acceleration (in the x-y plane) balances the component of 
the gravitational acceleration that is anti-parallel to it, leaving only an unbalanced 
^(-component of the gravitational acceleration. For a stationary axisymmetric space- 
time that has a reflection symmetry about an equatorial plane, in that plane there is 
no other component of the gravitational acceleration, so in the equatorial plane the 
local minima actually have zero acceleration and correspond to stationary geodesics 
or circular Keplerian orbits at the corresponding 
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For a nonstatic (e.g., rotating) stationary axisymmetric spacetime, the value hlo 
of the angular velocity that gives the local maximum of the acceleration (n_ < hlo < 
r2+) gives a local dehnition of a congruence, the SCAM, that in a local sense can be 
considered to be the most nearly at rest. Any slightly different rotation rate hi would 
give a change in the centrifugal acceleration and/or gravitational acceleration that 
would reduce the total acceleration. (In the Newtonian limit in which the gravita¬ 
tional acceleration is independent of the velocity, the reduction of the acceleration 
needed to balance gravity would be provided purely by the centrifugal acceleration, 
and in that limit, hlo = 0 is the angular velocity giving no centrifugal acceleration.) 

As one enters regions of strong gravity (e.g., near a black hole, r < 3M for a 
Schwarzschild black hole), one of the roots or 17+ may reach the corresponding 
endpoint hl-i or hl+i, or it may merge with hlo and thence go complex, in either case 
disappearing from the allowed region of the real 17 line (or all three roots may merge 
simultaneously in the nonrotating case), leaving only a single physical extremum 
(a minimum) for the acceleration. The fact that the acceleration then increases 
for a change in the velocity can be attributed to a reversal of the direction of the 
centrifugal acceleration |]^ |^, |^, |^, |^, |^, if one assumes 

(or defines) the gravitational acceleration to be independent of the velocity, or it 
can be attributed to a greater increase in the gravitational acceleration than that 
in the centrifugal acceleration [^, ^ 


@ 1 ^, 1 ^, 1 ^, 0 if one assumes 
(or dehnes) the gravitational acceleration to increase with velocity in the way (i.e., 
proportional to the square of the relativistic gamma factor, oc 7 ^ = 1/(1 — n^)) that 
one would get for an object moving transversely across a spatially flat horizontal 
floor in a rocket whose vertical acceleration in flat spacetime simulates gravity inside 
by the equivalence principle. 


3 Local definitions of stationary congruences 

The dehnition of the SCAM has the advantage of locality over the dehnition of the 
congruence with 17 = 0 (nonrotating with respect to inhnity), which requires the 
dehnition of which linear combination (ignoring the overall normalization, which is 
irrelevant for the present purpose) of the two Killing vector helds is k (usually made 
by choosing the combination that remains timelike at spatial inhnity), a dehnition 
that cannot be made locally but instead requires a knowledge of the behavior of the 
Killing vector helds out to spatial inhnity. 

The SCAM also has this same advantage, though to a lesser degree, over the 
ZAMOs (Zero Angular Momentum Observers), which are dehned to be orthogonal 
to the 1 Killing vector held and so require that that vector be uniquely picked out, 
again only up to normalization (typically by choosing the Killing vector held with 
closed orbits, which usually vanishes on a symmetry axis). This again requires 













nonlocal knowledge, unless one is at the symmetry axis where 1 vanishes. Since 
property (1) implies that 1 not vanish, I am explicitly assuming that one is not at 
a symmetry axis, except for some discussions below where I take the limit of going 
there. 

ZAMOs were originally called “locally nonrotating observers” because they 

have angular velocities midway between r 2 _i and ff+i, so that if two photons (in 
null but generically nongeodesic stationary orbits, say skimming along mirrors) were 
sent around both opposite directions from a ZAMO, they would both return to the 
ZAMO at the same time. This dehnition (essentially equivalent to dehning 1 to be 
the combination of the Killing vector helds with closed orbits) is quasilocal in that 
it does not require a knowledge of the Killing vector helds out to spatial inhnity, but 
it is still nonlocal in that it requires a knowledge of the helds along the stationary 
null orbits until they return to the ZAMO (i.e., all the way around the symmetry 
axis). On the other hand, the dehnition of the SCAM (when it exists) is completely 
local. 

To state more precisely what the conditions are for a quantity to be local if it 
depends on the two Killing vector helds k and 1, note that locally (away from spatial 
inhnity and from a symmetry axis) one has nothing that determines which vector 
in the entire k A 1 plane of vectors is k and which other vector is 1. In fact, one can 
make a global redehnition by the constant linear transformation 




(36) 


with constants a, /?, 7 , 6. 

If the four-velocity u given in Eqs. (pTf) , (pKl), and (^Tj) is to remain invari¬ 
ant under this redehnition of the Killing vector helds, the angular velocity hi must 
transform by the fractional linear transformation 


n 


n 


aVt — (3 
5 — yhl 


(37) 


Conversely, any four-velocity dehned in terms of an angular velocity hi which 
does not transform by Eq. (|T7|) , when the Killing vectors are transformed by Eq. 
(H, is not locally determined. For example, the congruence that is nonrotating 
relative to inhnity has hi = 0 when k is chosen to be the Killing vector held that 
remains timelike at spatial inhnity. But local information does not determine this k, 
and when one considers the transformation given by Eq. (|36|) , hi = 0 is not invariant 
under Eq. (^) as it would need to be to be a locally determined condition. 

The ZAMO is dehned by the condition that 


0 = u-1 = e-'*’(k -f Ozl)-l = e-'*’(-5 - O^C), (38) 
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or 


Qz = 


-6 _ _ 9t<fi 

C 


(39) 


ffii Qtptp 

However, this angular velocity, like hi = 0, does not transform according to Eq. 
when 1 is transformed to a linear combination of itself and of k, so it also does not 
locally determine a unique four-velocity u. 

To see this in detail and to prepare the way for recognizing which four-velocities u 
determined by the behavior of H, H, and C are locally determined and independent 
of the transformation (^61), one should note that Eqs. (p^) - (^) defining these 
quantities imply that they transform as 


A B 
B B 


a 


A B\ _/ a id\ / A B 
B c)~\j 5)\B Bj\f5 

/ a^A -|- 2af3B + j3‘^C a'-^A + {a6 + j3'y)B j35C 

VayH -|- (ah -|- (3^)B + (36C -|- 2'y6B -|- 


.(40) 


Therefore, if one makes a linear transformation of the Killing vector fields by Eq. 
(|36[) and then defines Clz = —B/C, one finds that when 7 7 ^ 0, this flz does not 
generally agree with the Qz one would get from applying Eq. (|37|) to Qz = —B/C. 

In contrast, other quantities and conditions are invariant under the transforma¬ 
tions of Eq. (^) and (^), such as u (by construction, since the transformation of 
n was designed to keep u invariant), a given by Eq. (|23|) , the normalized two-form 
(—obtained from the two-form A = k A 1 (which itself, like 
from Eq. ([T5|), gets multiplied by the determinant of the transformation matrix, 
ah — /hy), and the condition given by Eq. (p4D for the SCALEs. Thus the SCALEs 


are indeed locally determined congruences, unlike the nonrotating (relative to infin¬ 
ity) observers and the ZAMOs. 

Incidentally, I might point out that there are other quantities, analogous to A 
and that are not quite invariant under the transformations (|36|) and (pTl) but 
which for constant hi get multiplied by the constant (ah — /3y)/(h — yhl), namely K 
and e'*’ = F^/^. Thus $ = (1/2) In (-K-K) is not quite locally determined, but for a 
Killing vector field K that is fixed up to normalization (e.g., by being set parallel at 
some location to a stationary four-velocity u that is locally determined), V<h (which 
is the acceleration of u) is locally determined. 


4 Nonrotating Acceleration Worldlines (NAW) 
form a SCALE 

Another way of locally selecting a preferred stationary congruence of observers is to 
require that their acceleration vectors be Eermi-Walker transported and hence be 
nonrotating relative to an ideal system of gyroscopes carried by the observers. In 
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other words, the direction that each nonweightless observer feels to be np stays hxed 
in his or her Fermi-Walker transported frame, so that one might state Delphically 
that each such observer is “hxed up” or “stays up.” Such a congruence might be 
called Nonrotating Acceleration Worldlines (NAW). Here I shall show that a NAW 
is equivalent to a SCALE (or SCAM in the case that the local extremum of the 
acceleration is a local maximum). 

This equivalence is what one would expect in the Newtonian limit, since there, 
as discussed above, the SCAM has zero angular velocity (Hq = 0) and hence has 
an acceleration vector that stays constant. For observers that rotate, generically 
the acceleration vector also rotates, but for the SCALEs that are the two minima 
of the acceleration at angular velocities the acceleration vector is purely in the 
z-direction, which does not rotate. 

The Fermi-Walker derivative of a vector v along a worldline with four-velocity 
\i = d/dr and acceleration a = Vuu is [|| 

v' = Fu[v] = Vuv — (u A a)-v = Dw/dr — u(a-v) -|- a(u-v), (41) 

or, in component form, 

v'^ = U^Va-l3 - Uaa^V0 + (42) 

This is zero by construction when v is the four-velocity u. When v is the acceleration 
a, the Fermi-Walker derivative is 

a' = Fu[a] = Vua — a^u (43) 

(since u-a = 0). A Nonrotating Acceleration Worldline or NAW (more strictly, 
in view of my dehnition above of a NAW, a member of a NAW congruence) is a 
worldline in which a' is either zero or is parallel to a (i.e., a A a' = 0), so that in a 
frame carried along with the observer by Fermi-Walker transport, the acceleration 
vector a does not change direction. For a stationary observer, the magnitude of a 
stays hxed, so a' = 0 if the stationary observer is a member of a NAW. 

Using the explicit formulas ( p!7D and (^2]) for u and a with hxed H, which imply 

M(a;/3) + M(aa/3) = 0, (44) 

one can write this Fermi-Walker derivative of the acceleration of a stationary ob¬ 
server as 

a' = Vau (45) 

or 

= a^u^,0 = e-^^KPK^Kr^\KpK^,s - K^K0,s). (46) 
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( 47 ) 


Now a comparison with Eq. (|5B|) shows that 

oil 

Thus a' = 0 implies that the acceleration is an extremum with respect to hi, 
so NAW ^ SCALE (stationary worldlines whose acceleration is nonrotating also 
have their acceleration a local extremum). The proof of the implication in this 
direction used only the assumptions that the four-velocity is a linear combination 
e-'*’K = e-'*’(k -|- ni) of two Killing vector helds k and 1 that commute. Semerak 
[ 0 , after announcing that he had heard of my result by personal communication, 
has independently given a concise proof in this direction, NAW ^ SCALE. 

A mnemonic for remembering this result, using the somewhat cryptic shorthand 
phrases given above, is to say, “A stationary observer fixed up is extremely heavy,” 
or, “A stationary observer who stays up is extremely heavy.” 


5 A SCALE is a NAW 


To show the converse, that SCALE ^ NAW (stationary worldlines with locally 
extreme acceleration also have the acceleration nonrotating, which was not proved 
in [Q), we need to invoke not only the assumed commutativity of the two Killing 
vector fields k and 1, but also the assumed property (3), k A 1 A dk = kAlAdl = 0 
for the corresponding 1 -forms. 

For this purpose it is convenient to define the Killing vector field L (a linear 
combination of k and 1) that at the position of the stationary worldline under con¬ 
sideration has unit length and is orthogonal to the observer’s four-velocity u = e“^K 
there, 

L = yk -f 51, (48) 


where 7 and 5 are constants that at the position of the stationary observer have the 


values 


/9 




7 = - 


Vd 


G = 




-B-nc 

(B 2 _ AC){A + 2VLB + 112 ( 7 ) ’ 

A + VlB 

^(52 - AC){A + 2VLB + n'^C)' 


(49) 

(50) 


One can readily check from these formulas that at the position of the stationary 
observer, L-L = 1 and L-u = 0. Thus u and L are orthonormal vectors (or 1-forms, 
since I am using the same symbols for both vectors and the corresponding 1-forms, 
with the distinction, if necessary, being clear from the context) in the k A 1 plane. 
Also, clearly [k, 1] = 0 ^ [K, L] = 0, or 




( 51 ) 
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A 1-form in the (2,3) plane orthogonal to the k A 1 or u A L or (0,1) plane is 
a (which I will now assume is nonzero; otherwise it trivial that the acceleration is 
nonrotating). One may normalize it to get the unit 1-form 

a=(a-a)-^/2a (52) 

An orthonormal 1-form in this same (2,3) plane is 

b = *(u A L A a). (53) 

Thus (u, L, a, b) form an orthonormal frame or basis of 1-forms at the position 
of the stationary observer that is dehned purely in terms of the two commuting 
Killing vector helds and hence is Lie transported by the action of either one of 
them. However, this orthonormal frame is generically rotating relative to a frame 
that is Fermi-Walker transported along the stationary observer’s worldline. In fact, 
the Fermi-Walker derivative of any of the basis 1-forms above, or of any linear 
combination of them with constant coefficients, say v, is given by 

v' = *(u A a; A v), (54) 


where 

a; = - * (u A du) = * (K A dK) (55) 

is the normalized twist or rotation or vorticity 1-form of the Killing vector held K 
along one of whose integral curves the stationary observer moves. In particular, 

a' = *(u A a; A a). (56) 


The normalization of u makes it independent of the constant linear transforma¬ 


tion (|36D of the Killing vector helds and corresponding fractional linear transfor¬ 
mation (|37| ) of the angular velocity so long as the direction of K is kept hxed, 
as it indeed will be if one uses both transformations (^) and (|37|) . In particular, 
u is constructed to be independent of the normalization of K. Therefore, if the 
four-velocity u is locally determined, so is u and hence also a'. 

Eq. (^4|) is the analogue for 1-forms of the four-dimensional form of the equation 


that in the subspace of the tangent space orthogonal to the four-velocity u, the 
Fermi-Walker derivative v' is the three-dimensional cross-product of the rotation 
vector u with v. In component notation. 


/ 3 'y S 

Va = -eap^SU'^UJ^V , 


(57) 


where 




CUq. ^ ’ 2^®/37'5^ 


(58) 
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The Fermi-Walker rotation is thus in the two-plane orthogonal to the uAo; plane, 


so both u and a; have zero Fermi-Walker derivative. Also, it is clear from Eq. (^) 
or (1^ that oj is orthogonal to u. 

Note that in the definition (|^ or (^) of u, I am implicitly assuming that the 
four-velocity u which is differentiated there has a fixed angular velocity and hence 
is proportional to a fixed combination K = k -|- fll of the two Killing vector fields, 
with held constant during the exterior or covariant differentiation. In other words, 
u is the rotation 1-form for a congruence that has four-velocity everywhere parallel 
to a single Killing vector field K and hence may be considered to be rigidly rotating, 
since its shear and expansion are zero. 

When one considers a congruence, such as a SCALE or a NAW, that has 
varying as a function of the two coordinates (a = 2 or 3), then a; as I am 
generally using it in this paper is not the rotation 1-form that one would get from 
inserting that u(x“) into Eq. (|55[) , but rather the 1-form that one gets at each x“ 
by instead using the auxiliary congruence with constant (and hence with four- 
velocities parallel to a single Killing vector field K) that is chosen to match that 
of the original congruence at that point. (Of course, for a stationary congruence O 
does not depend on the proper time along each worldline. Unless the context makes 
it clear that I am assuming otherwise, I actually make the stronger assumption that 
the congruences I am considering have O and other scalar quantities independent 
of both of the (0,1) coordinates, so their derivatives in both the k and 1 directions 
vanish.) 

The stationarity of the worldlines tangent to the Killing vector field K imply 
that any covariantly determined vector v that is orthogonal to the worldline and 
determined by the local properties of the worldline and curvature (i.e., contractions 
of derivatives of the four-velocity and of the curvature tensor and its derivatives) 
is Lie transported by the Killing vector field K and so has itself a Fermi-Walker 
derivative that obeys Eq. (|5^) or (l57|). In other words, the stationarity of metric 
and of the worldline under translations along the integral curves of K imply that 
all locally determined vectors have components that stay constant in the locally 
rotating frame with basis 1-forms obeying Eq. ([^ ) or ( ^71) with respect to a locally 
nonrotating (Fermi-Walker transported) frame. 

If cu = 0, then one would obviously get a' = 0 by inserting a instead of v in Eq. 
so to show that SCALE 


(|54D or (1^ , so to show that SCALE ^ NAW, it remains only to check the case in 
which a; 7 ^ 0, which will henceforth be assumed. 

Now K A L oc k A 1, so Eq. (pT]) implies the corresponding equation for K and 

L, 


KALAdK = KALAdL = 0. 


(59) 


In terms of the rotation or normalized twist 1-form, the first of these equations is 
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simply 


a;-L = 0. 


(60) 


Therefore, u, a;, L, and 

L' = *(uAa;AL) (61) 

form an orthogonal set of 1-forms at the position of the stationary observer, with u 
and L having nnit magnitnde and hence with lo and L' having the same magnitnde. 
Thns one may write 

o; = =t:(u A L'A L), (62) 

and then nsing this in Eq. (p6|) gives 

a'= *(u A a; A a) = —(a-L')L. (63) 


iNOW mserrmg 


with Eq. (^H) shows that 


dsi _ {A + 2QB + Q‘^C) da 
27 dM ~ - AC 


(64) 


Pntting this expression into Eq. (^) and comparing with Eq. (^), or comparing 
with Eq. (^71) , one gets 


da^ ^ {A + 2QB + Q^C) da^ ^ 

4/D dn ~ 4VS2 - AC 'M 


(65) 


Thns, for a spacetime with two independent Killing vector fields with properties 
(l)-(3) above, da^/dfl = 0 a' = 0. That is, a Stationary Congrnence Accel¬ 
eration Locally Extremely (SCALE) is eqnivalent to a Nonrotating Acceleration 
Worldline (NAW). This proof implies the title statement, “Maximal Acceleration Is 
Nonrotating” (so it might be called the “MAIN Page proof.”) 

In contrast, one might be tempted to say of the converse resnlt given earlier, 
”NAW, it’s a SCAM.” Bnt this wonld be a scam, since actnally a NAW was proved 
to be a SCALE, and not all SCALEs are SCAMs. 

A Delphic way of stating the result, using phrases explained above, and substitut¬ 
ing “permanent” and “always” for “stationary,” would be, “Someone permanently 
fixed up is extremely heavy, and someone always extremely heavy stays up.” 


6 Acceleration and rotation of observers 

6.1 Stationary congruence parallel to a Killing field 

As an aside, one can make an analogy between the acceleration and rotation vectors 
for a stationary observer corresponding to a particular timelike Killing vector field. 
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and the electric and magnetic fields of an electromagnetic field as seen by a partic- 
nlar observer. In particnlar, if the timelike Killing vector field is K with sqnared 
magnitnde 

F = e^^ = -K-K, (66) 

then one may regard the normalized 2 -form 

f = —e“^dK = —du H—u A a (67) 

2 2 2 ^ ^ 

with antisymmetric tensor components 

fa/B — F (68) 

as being algebraically analogons to an electromagnetic field 2 -form with tensor com¬ 
ponents thongh it will not in general obey the analognes df = 0 and d * f = 0 
of the vacnnm Maxwell eqnations. 

(The normalization of f is chosen to make it independent of the constant trans¬ 
formations (|^) and (0) of the Killing vector fields k and 1 and of that keep K 
pointing in the same spacetime direction bnt which may change its overall magni¬ 
tnde.) 

Then in the frame of the stationary observer with normalized fonr-velocity u = 
the acceleration 1 -form with components 

= fa0U^ (69) 

is analogons to the electric field in a frame with fonr-velocity u, which has compo¬ 
nents 

Ea = FafjU^. (70) 

Similarly, the rotation 1-form with components 

UJa = (71) 

is analogons to the negative of the magnetic field in the observer’s frame, which has 
components 

5a = (72) 

One can then write the Fermi-Walker derivative of the components of any 1-form 
V which is Lie transported by the Killing vector field K as 

V'a = (73) 
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where 


= ( 5 / + + U''U0) 

= {5X + u^u^)f,,{5'^^ + u'^up) 

^o;/3 “ 1 “ 

fap ^a®/3 "h OjQ^Up ('^4) 


are the components of the rotation 2-form 

a) = — * (u A a;) = — ■^e'*’d(e“^'*’K) = —^e‘*’d(e“^u) (75) 

which is the part of the 2-form — ^du, or of the 2-form f = — ie“^dK, that is orthog¬ 
onal to u = e“'*’K. As usual, the normalization is chosen to make a) independent of 
the transformations (|3^) and (0- 

In particular, one can see from Eq. (^) that the Fermi-Walker derivative of 
the acceleration 1-form, a', is analogous to dvr times the Poynting vector flux of 
the electromagnetic held that is analogous to f. In other words, if one dehnes a 
symmetric second rank tensor 

47rt„/3 = (76) 

that is the analogue of the electromagnetic stress-energy tensor (generically not 
conserved here, since f does not obey the vacuum Maxwell equations), then 

a'„ = dvrpa = [*(uAa;Aa)]„ = -\-UaU^)t 0 ^v? = { 6 j-\-Uau'^)fp^f^^us. (77) 


= a-a 




One can easily see that there are precisely three independent Lorentz-invariant 
scalars at a point that depend only on the Killing vector held K and its hrst covariant 
derivative and are invariant under the transformations (|36D and (|T?|) that change its 
normalization: 

(78) 

(79) 

(80) 


U =U-U = /V 


and 


a-u- = = iiln/t A/t - h/“«/t)i‘''"- 


For example, dot products of a;, a and of all of its Fermi-Walker derivatives 
aO) of order n may be expressed algebraically in terms of these three scalars. For 
positive integers m and n. 


a'^ = a'-a' = = (a-a)(ci;-cu) — (a-o;)^ = a^o;^ — (a-o;)^, (81) 
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and 


= 0, 

(82) 


(83) 


(84) 

aO) .aO+2™-i) = 

(85) 


One can also get analogons relations between the vectors or l-forms themselves: 
If one dehnes 

a_L = a — (a-a;)a;/a;^, (87) 

the part of the acceleration vector or 1-form a that is orthogonal to the rotation 
vector or 1-form ut, which has sqnared length 

a]_ = a_L-a^ = o!'^ juP' = P — /oo'^, ( 88 ) 

then u, uj, a^, and a' form an orthogonal set of vectors or 1-forms, and each Fermi- 
Walker derivative a^"") is parallel (or anti-parallel, in alternate cases) either to a^ or 
to a': 


a{2n) _ 

(89) 

a(2n+l) _ (_l)-^;2r4a'. 

(90) 

When a' ^ 0, one can solve for the rotation 1-form 

oo = *(u A a' A a")/a'^ 

(91) 


in terms of the fonr-velocity u and the hrst two Fermi-Walker derivatives a' and a"of 
the acceleration a. In particnlar, the sqnared magnitnde of this rotation 1-form is 
then 


00 


2 _ 


= 00-00 = 


a" ■ a" 


a' a' 

[{AVB - BVA) + Vt{AVC - CVA) + 9?{BVC - 


cyB)f 


(92) 


One can compnte directly from Eq. (0) that this formnla for oo"^ in terms of A, B, 
and C is invariant under the transformations (j^) and (^T]). 

When one has a stationary observer in a curved spacetime, one can also get 
Lorentz-invariant scalars from contracting combinations of the Killing vector held 
K and its covariant derivatives with combinations of the Riemann curvature tensor 
and its covariant derivatives. As is well known, one can use Killing’s equation, Eq. 
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(^, to eliminate all covariant derivatives of the Killing vector field of order higher 
than two. Therefore, since 


f = —^e '^’dK = uAa + a) = uAa — *(u A a;), 


(93) 


any of these Lorentz-invariant scalars can be obtained by contracting appropriate 
combinations of u, a;, and a with appropriate combinations of the Riemann curvature 
tensor and its covariant derivatives, up to powers of the squared magnitude = 
—K-K of the Killing vector held that get multiplied by constants when one performs 
the transformations (|3^) and (^Tl), 

Just as a SCALE extremizes, and a SCAM locally maximizes, the scalar as a 
function of the angular velocity hi, so one could also dehne other congruences that 
extremize, maximize, minimize, or set to zero other combinations of the scalars a^, 

or a-a;. For example, one might dehne a Stationary Congruence Rotating At 
Minimum (SCRAM) as a stationary congruence that at each point, as a function 
of the angular velocity hi, minimizes the squared magnitude of the rotation 
u. (Remember that u is actually the rotation at each point, not of the original 
congruence being considered, but of an auxiliary rigid congruence, parallel to a 
single Killing vector held K with having a constant value that matches the x°'- 
dependent hi of the original congruence at the point where oj is being evaluated.) 
One can see from extremizing Eq. (|^) that the equation for this minimum, like that 
for a SCAM, is generically a messy quartic equation in O. The resulting four-velocity 
is independent of the transformations (|3^) and (l37|) and hence is determined purely 
locally. 


6.2 Nonstationary observers 

Many of the formulas above apply or have generalizations to the case of a nonsta¬ 
tionary congruence with u dehned over the region of spacetime under consideration, 
even when it is not proportional to any Killing vector held. For example, the ac¬ 
celeration vector a is still given by Eq. ([T^), and the rotation 1-form uj by the hrst 
expression on the right hand side of Eq. (|5^) . Then one can use the last expression 
on the right hand side of Eq. (p3| ) to dehne a 2-form f that is algebraically analo¬ 
gous to an electromagnetic held and which gives back the acceleration and rotation 
1-forms by Eqs. (|5^ and O- One can also use various other dehnitions, such as 
Eq. (|87|) for a^, the part of the acceleration vector that is orthogonal to oj. 

For a nonstationary congruence, the components of quantities dehned in terms 
of the four-velocity and its covariant derivatives (e.g., the components of a and of a;) 
are no longer constant in a frame which accelerates and rotates with the congruence, 
as such components are for a stationary congruence. Therefore, the Fermi-Walker 


derivative is no longer given by Eq. (|54D . In particular, a' is no longer given by Eq. 
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(|5B[ ) or (O)- Nevertheless, one can still define an analogue of dvr times the Poynting 
flux, 


dvrp = *(u A a; A a). 


(94) 


whose components are given by the right hand side of Eq. (fTTl) . When dvrp is 
nonzero, it may be normalized to define the unit 1-form 


L = p/(p-p)^^^ = [a^o;^ — (a-o;)^] * (u A a; A a). 


(95) 


Then u, a;/|a;|, a_L/|a_L|, and L form an orthonormal basis defined by the congruence. 

In general for a nonstationary congruence (or even for a stationary congruence 
in a metric which does not have a second Killing vector field obeying properties 
(2) and (3) above), L will not be parallel to any Killing vector. However, one can 
start from an original congruence with the four-velocity field u, construct the unit 
orthogonal vector L by the procedure above (when dvrp ^ 0), and thereby define a 
new set of congruences with four-velocities 


U 


u -I- uL 


(96) 


where v is the velocity of the new congruence relative to the old one, a parameter 
analogous to H in the stationary axisymmetric case. 

Then one can play the same game of choosing the velocity v at each point of the 
region of spacetime under consideration to extremize, maximize, minimize, or set 
to zero various combinations of such scalars as a^, o;^, a-o;, a-a', a;-a', a'-a', a-a", 
cu-a", a'-a", a"-a", etc., an arbitrarily large number of which can be independent for 
a nonstationary congruence, unlike the case of a stationary congruence, for which 
these scalars are all algebraically dependent on the first three. 

For example, one might define v to give a local extremum of a^, which might 
be called a Congruence Accelerating Locally Extremely (GALE), or a Congruence 
Accelerating Maximally (CAM) if the extremum is a local maximum. Or, one might 
minimize o;^, giving a Congruence Rotating At Minimum (CRAM). One could still 
have a Nonrotating Acceleration Worldline (NAW) if a' is parallel to a, but since this 
requirement is equivalent to solving two equations, it generically cannot be satisfied 
by any choice of the single parameter v for a nonstationary congruence, as it could 
in the stationary axisymmetric case. One could instead choose v at each point to 
minimize the squared length of the Fermi-Walker derivative of the direction of the 
acceleration, 

a , a , a^a'^ — (a-a')^ 


giving a Congruence Having Acceleration Rotating Minimally (CHARM). Clearly a 
NAW is a special case of a CHARM. 
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Because of the independence of the higher covariant derivatives of the four- 
velocity of a nonstationary congruence, in contrast to the case for a stationary 
congruence, it appears that one would not in general be able to conclude that a 
CAM is a CHARM, for instance. It might be of interest to see when this is so, but 
I shall not pursue this issue further here. 

Some of the formulas above can be applied even when there is no congruence 
but only a single worldline. In this case one cannot dehne a rotation 1-form lo by 
Eq. (^), but one can still dehne the acceleration a and its various Fermi-Walker 
derivatives a', a", etc. If one liked, one could make up a dehnition for a rotation 
1-form O', such as Eq. (HD, which agrees with Eq. (^5l) for a stationary congruence 
when 7 ^ 0. It might be somewhat better for a general stationary congruence to 
replace a' by its part that is orthogonal to a, namely 

a'j_ = a'— a(a-a')/a^, (98) 

with squared magnitude 

(a'_i_)^ = a'^-a'^ = a'-a' — (a-a')^/a^, (99) 

thus leaving out a possible change in the magnitude of the acceleration that is 
automatically absent for a stationary worldline, and to replace a" by (a'^)', the 
Fermi-Walker derivative of a'^. Then one could dehne a rotation 1-form to be 

a; = *(u A a'^ A (a'_L)')/(a'_i_)^ (100) 

One might prefer also to take only the part of (a'^)' that is orthogonal to a'^, but 
this makes no diherence in Eq. (|10(]| ). 

Thus for a single worldline one can say that it has the property of being a member 
of a NAW if and only if a'^ = 0, but without having the rest of a congruence, one 
cannot say whether or not it is a member of a SCALE, SCAM, SCRAM, CALE, 
CAM, or CHARM (except that it must have the latter property if it is a NAW). 

7 Stationary Congruence Accelerating 
Maximally (SCAM) 

For a generic location (of the two coordinates parametrizing the two-surfaces 
orthogonal to the two Killing vector helds k and 1) in a generic axisymmetric met¬ 
ric, the quartic Eq. (|3lD giving the angular velocity flo of the SCAM (Stationary 
Congruence Accelerating Maximally) is messy to solve, and the explicit form of the 
solution is presumably not very perspicuous (though I have not bothered to write 
it out in gory detail). However, there are at least three situations in which one can 
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get simpler explicit solutions. Below I shall give an example where the the quartic 
factorizes into two fairly simple quadratics. But hrst I shall give the angular velocity 
of a SCAM in a spacetime with relatively slow rotation. 


7.1 In a spacetime with slow rotation or weak fields 


For a spacetime with relatively slow rotation, after a suitable constant linear trans¬ 
formations (^61) of the Killing vector helds k and 1 if necessary, the part of the 
metric in the (0,1) plane is nearly diagonal, and so are its derivatives. In particular, 
-C —AC, and VB-VB -C —VA-VC. Then one can readily write down the solu¬ 


tion of Eq. (l34|) that is hrst order in B and VH, ignoring higher order corrections 
in these quantities (which will start with terms cubic in these quantities, which I 
shall label simply as 0{B^)): 

VA-{AVB - BVA) 


fin — 


0{B^ 


( 101 ) 


VA-(/lVC-CV/l) 

In the far-held limit outside an isolated source, typically A = —gu is just slightly 
less than 1 (e.g., roughly 1 — 2M/r at r S> M) and has a gradient nearly in the 
positive radial direction with coordinate r = x‘^. For simplicity, in the following 
discussion let the derivative with respect to r be denoted by a prime. (The context 
will keep it clear that here this prime does not mean a Fermi-Walker derivative, which 
will usually be zero for the quantities that I shall be applying the prime to, such as 
the stationary scalar quantities A, B, and C.) C = —g^pip will be negative, typically 
going roughly as the negative square of the distance from the axis of rotation, so at 
angle 9 = from the positive axis, C will go like —sin^ 9, and hence off the axes 
will have C < 0. Furthermore, AC' — CA' will be dominated by its hrst term and 
hence be negative. Then if the source is rotating in the positive (p direction, then in 
the far-held limit B = —gtip will be positive, typically going as (2 J/r) sin^ 6*, where 
J is the intrinsic angular momentum of the source Q. Thus B will have a negative 
radial gradient B' as it tends to zero at inhnity, and AB' — BA' will be dominated 
by its hrst term and hence be negative. 

Therefore, in the far-held limit outside an isolated source centered at the origin 


of standard spherical polar coordinates, Eq. ( |101|) gives the approximate angular 
velocity of the SCAM as 

B' 

This is in contrast to a ZAMO, which has 

That is, a ZAMO corotates with the source by what is usually called the dragging 
of inertial frames, but a SCAM counterrotates. Why is this? 


J 

^3 

(102) 

2.7 

■ 

(103) 
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A simple physical argument that applies in the far-held limit is the following: 
Consider the source being a ring of nonrelativistic matter in the the equatorial plane 
rotating with n 1 in the positive 9 ?-direction, and consider a stationary observer 
just outside this ring. In the observer’s frame, the energy density of the nearby 
portion of the ring will be a minimum if the observer is corotating with the ring. If 
the observer is not rotating relative to inhnity (fl = 0 ), she will see a higher energy 
density for the source, because of a relativistic 7 -factor (7 = 1/\/l — for the 
energy of each particle of the source (increase of its kinetic energy), and because 
of another 7 -factor from the Lorentz-contraction of the source. If the observer is 
counterrotating, she will have an even greater velocity relative to the nearby portion 
of the ring, and so she will see an even higher energy density (increasing roughly 
linearly with —when it is small compared with the magnitude of the rotation rate 
of the ring itself). 

Because energy density is the main source of the gravitational held in the nonrel¬ 
ativistic limit, a counterrotating stationary observer will thus have a greater grav¬ 
itational attraction to the ring (greater acceleration radially outward as seen in a 
freely falling frame). 

As the observer increases her counterrotation rate, she will also have a counter¬ 
acting centripetal acceleration inward, but it will increase at hrst only quadratically 
with her counterrotating angular velocity —Thus a small counterrotation rate will 
increase the experienced outward acceleration against gravity more than the increase 
of the centripetal acceleration inward, so the net acceleration (which is outward, 
primarily against the gravitational attraction of the ring) will at hrst increase with 
increasing the magnitude of the counterrotation angular velocity —(De Felice 
previously noticed this ehect on the equatorial plane of the Kerr metric.) 

Eventually, the (roughly quadratic) increase in the inward centripetal accelera¬ 
tion will balance the (roughly linear) increase in the outward acceleration against 
gravity, and one will reach a local maximum of the acceleration, thus at a counter¬ 
rotating angular velocity (flo < 0) of the SCAM. 

This is the picture in a weak gravitational held, but in a strong held (M ~ r), 
the increase of the gravitational attraction with velocity (by the two 7 -factors, i.e., 
as 1/(1 — v^)) can dominate over the increase in the centripetal acceleration (which 
goes as u^/(l — u^), but typically with a diherent coefficient, roughly 1 /r rather 
than roughly M/r^), so that there is no local maximum of the acceleration and 
hence no SCAM at that location. By taking out the 7 -factor dependence of the 
gravitational acceleration and including it instead with the centripetal acceleration 
(which is not strictly forbidden, since only the sum of the two is a gauge-invariant 
observable quantity), Abramowicz and his collaborators ||TI| , |T^, 0, n, 0,0, 

label this latter effect as the reversal of the sign of the centripetal 
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acceleration, though to me it seems more intuitively understandable to recognize it 


23 















as the greater increase of the gravitational acceleration with velocity than has the 

M, m m w 


centripetal acceleration 


31, 32 


7.2 At an orbital 

Another situation in which the SCAM has properties that are simpler than those 
given by the solution of the generic quartic Eq. is at a location which I shall 

call an orbital. This I dehne to be a location where there are two stationary time- 
like geodesics at different values of (e.g., corotating and counterrotating circular 
Keplerian orbits). There the angular velocity of the SCAM is given by a quadratic 
equation rather than a quartic. 

In terms of the (reversed-sign) {t,(p) metric components A = —gu, B = —gt^p■| 
and C = —g^pip dehned above, Eq. (^) gave 

VA + 

“ = 2(A + 2aB + av) ■ 

At a generic location, VA, Vi?, and VC are rather arbitrary vectors in the 2- 
dimensional (2,3) tangent plane, so for no value of is a = 0. 

Only if —2Vi? is on the two-branched hyperbola 0“^ VA-|-OVC is there a value 
of O that gives a = 0 (a stationary geodesic). If there are two stationary geodesics 
at different values of O, then the two branches of the hyperbola must degenerate to a 
single straight line (i.e., VA and VC must point in the same or opposite direction). 
For Vi? to be on this degenerate hyperbola (at two different values of O), it must 
also point in that same direction. In other words, at an orbital, VA, Vi?, and VC 
must all be parallel (or anti-parallel). However, this is a sufficient condition for an 
orbital only if both values of H that give a = 0 give timelike four-velocities. 

Since VA, Vi?, and VC lie in the tangent (2,3) plane that orthogonal to the 

plane of the two-form A = k A 1 (which is not to be confused with the scalar 

A = — k-k), the condition that these three gradients are all parallel is equivalent to 
the two scalar equations 

* (A A VA A VH) = 0, (105) 

*(AA VHA VC) = 0. (106) 

In the (2,3) coordinate space, these two equations generically have solutions only at 
isolated points, if at all, so it is by no means guaranteed that a stationary axisym- 
metric metric will have any orbitals. But if the space is symmetrical with respect to 
reflection about an equatorial plane containing 1, then on that plane VA, Vi?, and 
VC will all he within its tangent plane and be parallel (since they are all orthogonal 
to 1). Therefore, if a spacetime has such an equatorial plane, and if both values of 
H that give a = 0 give timelike four-velocities for some region of that plane, then 
this region consists of orbitals. 
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Assume that the common direction of VA, V-B, and VC is not orthogonal to the 
radial direction with coordinate r = x^. (he., choose the radial direction so that it is 
not orthogonal to that common direction.) As above, let the derivative with respect 
to r be denoted by a prime. Then clearly VA, Vi?, and VC will all be proportional 
to A', i?', and C respectively, with the same constant of proportionality (at a hxed 
location). 

Therefore, at an orbital, Eq. (|3^ for the angular velocity hi of a SCALE becomes 
[A' + 2B'Vl + C'fi^] [{AB' - BA') + (AC' - CA')fi + {BC - CB')n^] = 0. (107) 


This thus factorizes into two quadratic equations. One is the equation for the 
geodesic condition, a = 0, or 

A' + 2B'VL + C'O^ = 0. (108) 


Its solutions are the angular velocities of the Stationary Congruences Accelerating 
Locally Extremely (SCALEs) that have local minima of the magnitude of the ac¬ 
celeration (here both at the global minimum of 0). I have been calling the angular 
velocities of the SCALEs that are local minima but here I shall call them VL±Ki 
since for a stationary axisymmetric spacetime they are the angular velocities of 
circular Keplerian orbits: 




B' ± 


(109) 


The other quadratic equation. 


{AB' - BA') + {AC - CA')fl + {BC - CB')C^ = 0, (110) 


is the equation at an orbital for the Stationary Congruence Accelerating Maximally 
(SCAM), the SCALE with a local maximum of the magnitude of the acceleration 
(and for an unphysical root). The physical root has angular velocity 


fin — 


AC" - CA' + ^{AC - CA')"^ - 4{AB' - BA'){BC' - CB') 

2{CB' - BC) ' 


( 111 ) 


One can see from Eq. (|9^) that at an orbital, a SCAM has zero rotation, C = 0, 
so it is obviously also a SCRAM there. (This fact was previously noted for the 
Kerr metric, and now we see that it is general.) Therefore, at an orbital, a SCAM 
has the same four-velocity as a rigid congruence with four-velocity parallel to the 
Killing vector field K with constant O that at that location has zero vorticity. (The 


value of n giving a; = 0 on the equatorial plane of Kerr was first given in |2^, and 


later de Felice discovered this is also the angular velocity that extremizes the 
acceleration there.) 
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Again one should be reminded that I have dehned the rotation lo so that it is the 
rotation of a congruence moving along the orbits of a single Killing vector held K 
with constant (with the constant matching the angular velocity of the original 
congruence at that location, but not necessarily at other locations). That is, if one 
inserted the ^“-dependent four-velocity u(x“) of a SCAM at an orbital into Eq. (|55D , 
a; = ^ * (u A du), one would generically not get zero, but only if one used in that 
formula the four-velocity u = e“'*’(k -|- fll) with constant chosen to make this u 
match that of the SCAM at the position where lo is being evaluated. Thus a SCAM 
at an orbital is not generically part of a congruence that itself has zero vorticity. 
Only the rigid congruence moving along the orbits of K has zero vorticity there (and 
that congruence generically has nonzero vorticity at other locations). 

(When we have two Killing vector helds obeying the properties (1) - (3) above, 
stationary congruences with zero vorticity are those with a constant ratio of angular 
momentum to energy. 


angular momentum Uu 


energy 


—Uo 


-B-CVL 
A + BVL 


= const.. 


( 112 ) 


such as the ZAMOs with kl obeying Eq. ( p.03|) , kl = —BjC^ so that its angular 
momentum, and hence the ratio above, is zero. Thus a congruence of ZAMOs has 
the local property of zero vorticity |^, though it generically does not have 
a; = 0 by my indirect method of dehning lo in terms of an auxiliary congruence for 
each location that rotates rigidly with constant kl. Incidentally, though the zero 
vorticity of a ZAMO is a local property of that congruence, it is not sufficient to dis¬ 
tinguish locally that particular congruence from the other zero-vorticity stationary 
congruences with different constant values of the ratio (|112|) of angular momentum 
to energy. The problem is that locally there is no way to distinguish the angular 
momentum connected with the Killing vector held 1 with closed orbits from a com¬ 
bination of angular momentum and energy connected with a different Killing vector 
held, since the property only 1 has, of having closed orbits, is not a local property 
that can be determined without knowing the metric in a loop around the symmetry 


axis. 


At an orbital one has, for a hxed choice of the Killing vector helds k and 1, seven 
special values (at least) of the angular velocity kl: kl+i and k2_i given by Eq. (|35|) 
(those of the speed of light in the forward and backward directions respectively), 
fl+K and fl-K given by Eq. (|109|) (those of the stationary geodesics or circular 


Keplerian orbits in the forward and backward directions respectively), klo given by 
Eq. ( |111|) (that of the Stationary Congruence Accelerating Maximally, or SCAM, 
which has the local maximum of the magnitude of the acceleration as a function of 
kl, and which, like the geodesics with Q±k that are members of the other Station¬ 
ary Congruences Accelerating Locally Extremely, or SCALEs, are also Nonrotating 
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Acceleration Worldlines, or members of a NAW congruence that Fermi-Walker trans¬ 
ports the acceleration vector), Qz given by Eq. (that of the ZAMO, which has 
zero angular momentum), and ^2]VR = 0 (that of a congruence nonrotating relative 
to inhnity). 

The four-velocities corresponding to the hrst hve of these angular velocities are 
determined locally and are invariant under the linear transformations (^^ of the 
Killing vector helds k and 1, but, as discussed above, that is not true of the last two, 
since they require the nonlocal knowledge of which Killing vector held is k (e.g., 
the one that is timelike at inhnity, usually normalized to have unit timelike length 
there) and of which one is 1 (e.g., the one with closed orbits, usually normalized to 
give period 27r around the orbit). (For the hrst two angular velocities there 
are no normalized four-velocities with those angular velocities, so for them I mean 
instead the corresponding null vectors, which are locally determined only up to 
normalization.) Of course, the particular values of all but the last of these angular 
velocities O depends on the particular choice of k and 1, but my point is that the 
corresponding K (up to normalization) at each point does not, for the hrst hve 
angular velocities. 

Now, as one might expect, there are a number of algebraic relations between these 
angular velocities and between the corresponding four-velocities. For example, it is 
well known [Q that the ZAMO has the average of the angular velocities of the two 
null orbits, 

= 2(^+1 + ^-i)’ (113) 

which implies that if a ZAMO sent two photons around an orbital in opposite di¬ 
rection (using, say, a tube to dehect the photons into these nongeodesic null orbits 
by an inhnite number of glancing collisions that each transfer an inhnitesimal mo¬ 
mentum to the corresponding photon), they would return to the ZAMO at the same 
time [0. Of course, this procedure requires the nonlocal information of the metric 
around the orbital, the same nonlocal information that is required to pick out 1 as 
the Killing vector held with closed orbits, a choice that is necessary before one can 
dehne the angular velocities and get formulas such as Eq. (tll3 )i which is not invari¬ 
ant under the transformations (|36D and (^) if there the transformation constant 
7 7 ^ 0. However, one may at least note that Eq. (|113|) is independent of the choice 
of the Killing vector k (determined by the transformation constants a and /3) or of 
the normalization of the Killing vector 1 (determined by the transformation constant 


5 if 7 = 0), even though changing these will change the H’s appearing in Eq. ( 113 ). 
Another more complicated relation one may hnd is 


\/ (f^+i — ^-k)(^-k — + y (f^+1 — f^+x)(f^+x — H-i)H-ii' 


(H_i_i — — H_i) -|- W (H+i — H_|_x)(H+x — H-i) 
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Thus the SCAM angular velocity hlo at an orbital is a weighted average of the angular 
velocities ^±k of the two circular Keplerian orbits. This relationship is invariant 
under the transformations (^) and 

This relation simplihes greatly if one considers instead the relative velocities 
between these various observers. One readily hnds that in the frame of an observer 
having four-velocity ui with angular velocity Oi, the three-velocity of a second 
observer having four-velocity U 2 with angular velocity ^2 is 


V 2 = 


U 2 


-U 1 U 2 


- Ui = 


(115) 


where the (signed) relative speed of observer 2 relative to observer 1 is |34 


_ AC{yt2 - Oi) 

A + B{Vti + Vt2) + CVtiVt2 

(o+i — r2_i)(r22 ~ Oi) 

(0+1 + 0_i)(02 + Oi) - 20+i0_i - 2 O 1 O 2 ’ 


(116) 


which is invariant under the transformations (|5BD and (P7|). 

Then one can calculate that at an orbital, the two circular Keplerian orbits 
have equal and opposite speeds in the frame of the SCAM, the magnitude of which 
may be called the Keplerian orbital speed vk'- 


Vk 


^^(O+ii:, Oq) — —u( 0 _i(', 0 o) 

\l (0+1 — 0_/f)(0+/f — 0_i) — (0+1 — 0+j^)(0_x — 0_i) 
Y^(0+i — Q-k){^+k — 0_i) -|- (0+1 — 0+i^)(0_^ — 0_i) 

(\/(^+i ~ 0_ii')(0+i;s' — 0_i) — (0+1 — 0+ii')(0_ii' — 0_i))^ 

(0+1 — 0_i)(0+i^ — 0_x) 

(^2 - Acy - y/(S2 - Acy^ - 4(52 - AC){B'^ - A'C) 
2 ^( 52 -AC)(5'2-A'C") 

D'-^ 2^/DH 

2y/DH ~ D' + 


(117) 


where Eqs. (^) and (|2^) give D = B^ 
two-dimensional block of the metric). 


AC (minus the determinant of the first 


H = B'^- A'C', 


(118) 


and 


a = D'^-ADH = {B^ - AC)'^-A{B^ - AC){B'^ - A'C') 

= {AC' - CA')^ - 4{AB' - BA'){BC' - CB') 

= C2C"2(0+1 - 0+^)(0+i - 0_^)(0+i^ - 0_i)(0_i^ - 0_i). (119) 


























One may check from Eq. for the transformations of the quantities 5, and 
C that the hnal three expressions on the right hand side of Eq. (|117|) are indeed 
invariant under the transformations (^) and (|37|) and hence are locally dehned. 

Therefore, a SCAM at an orbital has a four-velocity uq that is a normalized 
average between the four-velocities u+x and u_x of the two circular Keplerian 
orbits: 

u+x + u_x _ u+x + u_x 


Uo = 


|U+x + u_x| 




( 120 ) 


2 u+x-u_x 

This means that if one has two equal-mass point particles moving along opposite 
stationary Keplerian orbits, and they collide to form a single particle in a totally 
inelastic collision, the velocity of this single particle immediately after the collision 
will be that of the SCAM at that location. This relation was discovered by Semarak 

Thus we 


Erst in the special case of the Kerr metric , and then later in general 
see that it is a feature at an orbital of any stationary axisymmetric metric invariant 
under reversing both t and ip. 

One can see that even if VA, VB, and VC are all parallel (or anti-parallel), the 
SCAM is not dehned as a real congruence when a dehned by Eq. (|119|) goes negative. 
It reaches zero when one of the two circular Keplerian orbits reaches the speed of 
light (usually when Q-k becomes as negative as O.i for a positively rotating source, 
so that it is usually the last factor in the last expression of Eq. (|119|) that goes to 
zero hrst). Then there are no longer two timelike stationary geodesics, so one is not 
really at an orbital as dehned above, even if VA, VB, and VC are all parallel (or 
anti-parallel). 

At an orbital, the two SCALEs that are not the SCAM give the two circular 
Keplerian orbits. Away from an orbital, these two SCALEs have locally minimal, 
but not zero, magnitude of acceleration. One might think that Eq. (|120|) would 
generalize to this case, with u±x being replaced by u±, the four-velocities of the 
SCALEs that are not the SCAM, but this is not generically the case. In other 
words, Eq. (|120|) applies only at an orbital, the only place where there are stationary 
geodesics, the circular Keplerian orbits with four-velocities parallel to combinations 
of the two Killing vector helds k and 1. (I am always implicitly excluding stationary 
geodesics that have four-velocities parallel to any other possible Killing vector helds 
that might be present, such as in a spherically symmetric spacetime.) 

Another diherence from the case of of an orbital, where the solution of the 
quadratic equation for the SCAM breaks down (becomes complex) after one of the 
circular Keplerian orbits reaches the speed of light as one enters deeper into a strong 
gravitational held, is that away from an orbital, the solution of the quartic equation 
for the SCAM becomes complex after the four-velocity for the SCAM merges with 
that of one of the other SCALEs as one enters deeper into a strong gravitational held. 
(Deeper in the held there is only one extremum, a minimum, for the magnitude of 
the acceleration as a function of the angular velocity for timelike worldlines.) Thus 
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one can have the case in which the relative velocity between the SCAM and one of 
the other SCALEs goes to zero (at the boundary of the region where the SCAM is 
dehned as a real congruence), whereas the relative velocity between the the SCAM 
and the third SCALE can remain nonzero there. 

In the frame of the SCAM at an orbital, one can readily show that the 3-velocity 
of an observer nonrotating relative to infinity (i.e, nonrotating relative to the k 
Killing vector, so hlATK = 0), and of a member of a ZAMO are, respectively. 


vnr — 


-noVB^ - AC 

A + Bilo ’ 


B + CQr 


Vz = 


( 121 ) 


( 122 ) 


VB’’- - AC 

However, these quantities are not locally dehned and hence are not invariant under 
the transformations (|36|) and (^T]). 

For a corotating source, usually A, B, —C, and —hlo are positive, so both vnr 
and Vz are then positive. In the far-held limit outside an isolated source centered 
at the origin of standard spherical polar coordinates, one gets, using Eqs. (|1U2|) and 


(HID, 


Vnr - 
BC 


^/^B' 
C' 
■CB' 


J • a 

— sm 0, 


Vz 


• ft 
— sm 6 


3vnr- 


(123) 

(124) 


The orbitals of an isolated source, if any exist, are on or near the approximate 


equatorial plane 9 = 0, but Eqs. (|123|) and ( p.24|) apply at arbitrary 9 outside an 
isolated source in the far-held limit. 

The stationary observer of the SCAM at an orbital has an acceleration that, in 
an orthogonal {g 23 = 0) coordinate system for the (2,3) plane in which the gradients 
of A, B, and C are purely in the direction of the coordinate r = x^, has magnitude 


Oq 


4:D 

{B^ - Acy - ,J{B^ - Acy^ - 4(^2 - AC){B'^ - A'C) 
A{B^-AC)^r 


(125) 


(Here the subscript 0 does not denote the time component of the acceleration, which 
is zero, but rather the acceleration at the SCAM, for which the subscript 0 has been 
used.) The acceleration 1-form of the SCAM is then a = , with 


e'’ = = y/^Vr 


( 126 ) 
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being the unit l-form in the outward radial direction. This acceleration l-form a is 
clearly locally dehned and invariant under the transformations (|36|) and (^). 

Now one can readily calculate that for any stationary observer (one moving in 
the k A1 plane along the orbits of a hxed Killing vector held K) that has speed v in 
the frame of the SCAM at an orbital, the acceleration l-form is simply 


a — oq 


1 — 

1 — 


(127) 


In particular, the acceleration is a symmetric function of the velocity in the SCAM 
frame. 

One can regard ao/(l —as being the (radial outward) gravitational part of the 
acceleration, which is indeed proportional to 7 ^ = 1/(1 — u^) as one would get by the 
equivalence principle for an object moving horizontally across the hat hoor of a rocket 
that has constant acceleration perpendicular to the hoor. Then [ao/v]^) vy{i-v^) 
can be regarded as the centripetal part of the acceleration, pointing radial inward, 
and also having the 7 ^ dependence that it does in hat spacetime. 

For example, suppose that one had an idealized (hctitious) static, spherically 
symmetric metric with constant m. 


ds^ = sin^ Odcp^ + dR^ + R^dO^, (128) 


which is spatially hat but has a spherically symmetric gravitational potential 

$ = (129) 

for a static observer with four-velocity u = e“'*’K = e“'*’k. (I use the radial coor¬ 
dinate R instead of r so that later I can compare with a diherent coordinate r in a 
diherent metric, such as Kerr-Newman. The radial metric component is qrr = 1, 
so the unit radial l-form e'’ is simply dR or Vi?.) 

Here B = gtip = 0, so the SCAM consists of nonrotating (static) observers, which 
have acceleration 

TTl 

a = V<F = —e'’= ooeC (130) 

R^ 

In the equatorial plane, 9 = 7 r/ 2 , which has orbitals everywhere, a stationary ob¬ 
server orbiting with velocity v relative to that of the static SCAM has acceleration 


a — Un 


ac = ( 


m 


i?2(l 


i?(l 


)e' = ao- 


vVvl, 


(131) 


with oo = m/R^ and vk = \/aoR = \Jm/R, so it is natural to split up the total 
acceleration into a gravitational piece 


-e = 


Oq 


m 


a.g 


i?2(l 


(132) 
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and a centripetal piece 


a, = 


R{1 


-e = 


aov‘^/vK 


One might nse this analysis with the idealized metric 
orbital radins of cnrvatnre’ 


_ '^K V^) y/9rr 4/^^ g^r 


Oo 


and an ‘effective gravitational mass’ 


m = — = Rv\ = 
Oo 


H 


D' + y/^ 


(133) 

to dehne an ‘effective 

(134) 


R^H (D' - VF)i^ 

QrrD H 


(135) 


for an orbital in any stationary axisymmetric metric. These formnlas are chosen to 
make Eqs. (|131|) - ( p.33|) trne for a particnlar orbital in any stationary axisymmetric 
metric. Then one can write the total acceleration (|127]) as 


a = 




1 


e 


(136) 


which is thus much simpler in the SCAM frame than the corresponding formula in 
the ZAMO frame H- 


Following [^, we can also note that there is a simple geometrical description 
of the effective orbital radius of curvature R\ If the Killing vector held 1 has period 
A(p = 27r, then the circumferential radius (circumference divided by 2ti) measured 
by a stationary congruence with constant is 


= Jd/F = 


(137) 


Then if 

dSr = y/^dr 


(138) 


is an inhnitesimal element of proper distance in the radial direction (an element of 
proper radius), one can readily calculate that if = ffg (to be held constant during 
the spatial differentiation, as usual in this paper). 


dSr / \ 

R = f^, (139) 

which is what one would calculate the proper radial distance to be to a center (where 
f should vanish) if one assumed that f varied linearly with proper distance. E.g., for a 
circle of latitude in the northern hemisphere on the surface of an axisymmetric earth 
in hat space, f would be the cylindrical radial distance from the circle to the axis 
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inside the earth (intersecting it orthogonally), dsr would be an infinitesimal proper 
distance along a meridian (line of constant longitude) in the southward direction 
along the earth’s surface, and R would be the distance along a cone tangent to the 
surface of the earth at the circle, from the circle to the apex of the cone over the 
north pole. 

Then if one uses Eqs. ( |126|) , ( p.36|) , ( |138|) , ( |139|) , and the fact that df = Vf = 
{df/dsr)e'^, one gets an alternative simple formula for the acceleration of a stationary 
worldline of speed v in the frame of the SCAM, 


a = 


Jr 


V$ —n^Vlnr 




1 — / r 


(140) 


Of course, unlike the case for the idealized metric (|128|) , for a more general metric 
not only the effective orbital radius of curvature R, but also the effective gravita¬ 
tional mass m, may vary from orbital to orbital. For example, for the Schwarzschild 
metric 


ds^ = —(1 — ‘^^)dt^ -I- sin^ Odif"^ + (1 ~ 

r 

with constant ADM mass M, one has 

ao = ^(l-) 


2M, 


)-Mr^ + r^de^ 


(141) 


(142) 


and 


M , 
= \/—(1 


2M, 


- 1/2 _ 


M 


/y /y \I T' 

SO one gets for the effective radius R and mass m 

2M. 


2M’ 


R = r(l — 


r 


1 - 1/2 


and 


, 2M,_o/2 

m = M{1 -) 


(143) 


(144) 


(145) 


In the Schwarzschild metric, one can compensate for this dependence of the 
effective mass m on the radial coordinate r by defining an ‘effective Scharzschildean 
mass’ 


M = 


m 


Vr 


(1 + 2m/R)-^R ao(l + 2r;|)3/2 
[(D' +v^)2 + 8M]3/2’ 


(146) 


which for the Schwarzschild metric is designed to give precisely M. 

Since only the total acceleration a is observable, one is free to divide it up into 
‘gravitational’ and ‘centripetal’ contributions any way one wishes. Although the 
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split given by Eqs. (|131j) - ( |133|) above seems most natural to me (and to several 
others ^ ^ ^ ^ ^ ^ Abramowicz and his collaborators 

ini , [T^ , 1^, 14, ]^,|I6,|I7,^, ^ have advocated an alternative split 


cl — cl^ 3.^ 

(147) 

with 



(148) 

being independent of velocity and all the velocity dependence being put into the 

centripetal piece 

ao{l-v]^)v‘^ ^ 

2/1 2^ ^ ■ 

(149) 

If one wished, with this alternative split one could define 
radius R and effective gravitational mass rfi so that 

an alternative effective 

rfi 

ag = aoe^ = 

(150) 

and 

2 


V 

Sir = — =- 

i?(l -^2) 

(151) 


Then the alternative effective radius is 


R — ^ ^ ^ 

ao{l — v]^) 1 — vj^ R — m' 


and the alternative effective gravitational mass is 


m = 


V 


4 

K 


ao(l - 


mR^ 

{R-mf 


(153) 


One can also give the same geometrical interpretation for R as for R if one uses 
instead of f the optical circumferential radius 


from the optical metric 
Then one has 


and 


a = ooe'' 


1 — 7 - 



(154) 


(155) 

dsr 

dr 

(156) 

dh _ 

— = V<I> — -r I Vlnr. 

r \1 — J 

(157) 
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When the circular Keplerian orbital velocity vk exceeds 1 (the speed of light), 
Sc, dr, and R reverse sign, which is what Abramowicz calls the reversal of the 
sign of the centrifugal acceleration. However, this interpreted reversal depends on 
the particular splitting of the total acceleration a into gravitational and centripetal 
parts given by Eqs. (|147|) - (|149|) so that the gravitational part is dehned to be 
independent of the orbital velocity v. To me it seems more natural to use the split 
given by Eqs. (|131|) - (|133|) so that the gravitational part a^ has the 7 ^ velocity 
dependence that one would expect if gravity couples to energy rather than, say, to 
rest mass. In an equivalence-principle argument, one can readily calculate that the 
acceleration for objects skimming horizontally over a spatially flat floor of a rocket 
having constant acceleration vertically in flat spacetime indeed has this 7 ^ velocity 
dependence. 

Nevertheless, even in the rocket example, one could follow Abramowicz and 
say that part of the acceleration is centripetal acceleration from the curvature of 
the floor world membrane (the timelike three-surface that is the world history of 
the accelerating floor). Indeed, the floor looks curved as seen by photon geodesics 
(which the floor will intercept twice if they are not traveling purely vertically), even 
though at each moment of time the instantaneous two surface of the floor is flat 
and contains tachyon geodesics that move at inhnite speed. However, one cost of 
Abramowicz’s interpretation is that no gravitational acceleration would be ascribed 
to photons, so that they would be interpreted as having no weight, even though they 
have energy. 


7.3 Near a rotation axis 

A third location in which the properties of the SCAM are more simple than in the 
generic location (besides the two discussed above, hrst of the slow-rotation and/or 
far-held limit that one would expect to hnd far from an isolated source, and second at 
the location of orbitals where both corotating and counterrotating circular Keplerian 
geodesic orbits exist) is inhnitesimally near a rotation axis where 1 vanishes. If 1 
is given its standard normalization of having closed orbits with period 2n (i.e., if 
1 = d/d(p with ip being periodic with period 2n), then regularity of the metric near 
the axis implies that C = —IT must go as — -|- O(tu^) (with unit coefficient), to 

lowest order being the negative square of the distance vu from the axis along a spatial 
geodesic that is orthogonal to k and 1 and which intersects the axis orthogonally. 
We can let w be one of the cylindrical coordinates for the (2,3) plane orthogonal to 
k and 1 and let z be the other, with the property that z is proper distance along 
the axis and that each spatial geodesic mentioned above that intersects the axis 
orthogonally has constant z, so that Vcu and Vz are orthogonal everywhere and the 
former is normalized to have unit magnitude. 
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Then in a neighborhood of the axis where these two coordinates are well behaved, 
the metric may be written 


ds‘^ 


—A{zu, z)dt^ — 2B{w, z)dtd(p — C{w, z)d(p‘^ + dw"^ + gssizu, z)dz^ 

— [Aq{z) + A2{z)w‘^ + 0{w^)]dt^ — 2\p.z{z)w‘^ + 0{w'^)]dtdLp-\- [w‘^ + 0{w^)]d‘{)^ 
+dw‘^ + [l + 0{w‘^)]dz^. 

(158) 


On the axis the ZAMO is locally determined, since 1 is uniquely determined by the 
requirements of the previous paragraph if one has access to the axis, unlike the case 
away from the axis, where one does not have local access to the axis. In particular. 






= lim 


B{w, z) 
C(zu, z) 



(159) 


where the last expression is to be evaluated on the axis, at ci7 = 0. However, note 
that although the ZAMO is locally dehned at the axis, the angular velocity Qz 
ascribed to it is not, since the latter depends on the dehnition of k, which is not 
locally determined but can be redehned by multiplying it by a constant and by 
adding any constant multiple of 1 . 

Similarly, 

Ao(z) = A(zu = 0,z), (160) 

and 


Ay 


= lim 


VA(zu,z)-VC(zu,z) _ 4 

0 4C(w,z) 


2A''(z)-A'(zr/Ao(z)), (161) 


where in the last expression the (four-dimensional) Laplacian of A is to be evaluated 
on the axis, and I am here using the notation that a prime denotes a partial derivative 
with respect to Here I shall also use the convention that when I give the functional 
dependence of any quantity as (z), I mean that the quantity is to be evaluated at 
zu = 0. 

If = dz is the unit 1-form in the ^(-direction on the axis, then one can readily 
calculate that the ‘nonrotating’ (H = Htvj? = 0 ) congruence with four-velocity u = 
has, on the axis, an acceleration 1-form 


a(;2) = a{z)e^ 



and a normalized rotation 1-form 


^nr{z) = — —Aq 


(162) 


(163) 
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This is negatively rotating, since it is the ZAMO (and not, for example, the SCAM) 
which has the angular velocity of the Killing vector field K = k + hl^l, which, 
when VLz is held constant, has zero rotation or vorticity precisely on the axis, uiz{z) = 
0. (This is rather opposite to the case at an orbital, where it is the SCAM rather 
than the ZAMO which has the angular velocity of the Killing vector held that has 
zero rotation or vorticity there. Thus a SCRAM, which has O chosen at each location 
to minimize o;^ there, should interpolate between the SCAM at any orbitals and the 
ZAMO on an axis of symmetry.) 

_-j /o 

One should also note that precisely on the axis where 1 vanishes, u = (^)k 

for any O, so the acceleration a{z) there is independent of O. Thus, strictly speaking, 
the dehnition of a SCALE does not work precisely on the axis, but it does if one 
moves inhnitesimally off, and then one can dehne the SCALEs on the axis to have 
the limit of Vt{w) as ci7 —>■ 0. 

Now when one evaluates the quantity on the left hand side of Eq. (|^) for a 
SCALE near the axis of symmetry, one hnds that it goes as plus higher-order 
terms in as the acceleration goes as a(z) plus a correction term that to lowest 
order is proportional to with a coefficient that depends on hi. If one divides the 
left hand side of Eq. (|3^) by and takes the limit of w going to zero, one gets on 
the axis not a quartic but a cubic for the SCAM hlo and the other two SCALES 

(fi - nzf - (A2 + - ^)(f] - ^z) + \^^'z = 0. (164) 

Here all the quantities are to be evaluated on the axis zu = 0 and so are functions 
purely of z, and for compactness I have used A and its 2 ;-derivative A! instead of Aq 
and Aq, since they are the same quantities on the axis. 

There are three distinct real solutions of this equation, giving one local maximum 
(the angular velocity hlo of the SCAM) and two local minima (the angular velocities 
of the remaining two SCALEs) of the magnitude of the acceleration, if the 
discriminant of the cubic, which is proportional to 

A >2 

27A’^n’l - 16(A2 + - ^)^ (165) 

is negative. If, on the other hand, this quantity is positive, there is only one real 
solution, corresponding to the minimum value of the acceleration, and thus the 
SCAM does not exist there. 

One notes that the coefficient of the quadratic term in O — VLz is zero. Thus 
when all three solutions are real and hence denote the angular velocities of the three 
SCALES, 

Oo + + 11+ = 311z, (166) 

so the ZAMO angular velocity is the average of that of the three SCALES on the 
axis. 
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Instead of writing the cnbic Eq. (|164|) for the SCALEs in terms of their angnlar 
velocities (which depend on the choice of k), it may be more illuminating to write it 
in terms of the orthonormal 2 ;-component a;, on the axis, of the normalized rotation 
1-form uj of the rigidly rotating congruence with the corresponding angular velocity 
at the location in question on the axis, 

ljj{z) = uje^ = {z)(yt — Vtz)G^ ■ (167) 

Then a bit of algebra, including the use of the standard formula [Q 

'Ca;/37 = (168) 

for any Killing vector held ^ to eliminate its covariant derivatives of order higher 
than one, gives the following simple cubic for the orthonormal rotation component 
of the limit of a SCALE on the axis: 


u!^—puj — q = 0, (169) 

where 

P = RtU^ 

and ^ ^ 

Q = ■^<iR^vozi = 2«(^h + a^n). (171) 

Here all of the quantities (such as the orthonormal Riemann and Ricci curva¬ 
ture components, and the orthonormal 2 ;-component a of the acceleration—its only 
nonzero component) are to evaluated in the limit of going onto the axis, the prime 
denotes a derivative with respect to proper length z along the axis, and 

ujji = A — Vtz) ( 172 ) 


is the z-dependent rotation of any rigidly rotating (flij = const.) congruence, e.g., 
the nonrotating congruence with = 0, though to dehne this particular congruence 
requires a specihc choice of k that is not required in Eq. 


One can easily 

see that + clojr) is invariant under changing from one allowed rigidly rotating 
congruence to another by changing the constant in Eq. ( |172| ), so the solutions 
of the cubic Eq. (|169|) do not depend on this choice. In other words, the coefficients 
p and q are both invariant under the allowed transformations (|3^ ) and (|37|), which 
on the axis are restricted to have 7 = 0 and 5 = 1 so that only k, but not 1, may be 
changed. 

The explicit solutions of the cubic Eq. (|169|) may be written as 


o /P • fl • -1 

UJO = -2Y'-sm[-sm 


iQ\ 
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4p3 


(173) 











for the normalized rotation rate of the SCAM, and 


[p . r 271 1 . _i / / 27 

.,* = 2y/-sm[±---sm- (,y^)|. 


( 174 ) 


for the rotations of the other two SCALEs on the axis. 

After solving the cubic Eq. (|169|) for the invariant normalized rotation rates oj, 


one can use Eq. (|167| ) to solve for the angular velocities hlo and if one wishes, for 
a particular choice of k and hence of A and of Viz- However, these angular velocities 
n are not invariant under the allowed transformations (p6D and (|37|) that change 
k, whereas the normalized local rotation rates u (of congruences with constant 
angular velocities hi that match those of the SCALEs at the chosen location on or 
inhnitesimally near the axis) are invariant under these transformations of the Killing 
vector held k. 

To hrst order in g(27/4p^)^/^ when the explicit solutions (|173|) and 

(|174|) reduce to the approximations 


Wo ~ — = -B -^- 

P R 


t UU t'UJ 


— 


W± 




(175) 


(176) 


For example, in the far-held limit outside an isolated stationary source of mass 
M and intrinsic angular momentum J at rest at the origin, with z ^ M + ^fj 
being the positive proper distance along the axis from the source, in the direction 
of the angular momentum vector, and with the stress-energy tensor (and hence 
7?££, assuming Einstein’s equations) being negligible there, one has on the axis A = 
—gtt ~ 1 — 2M/z, a = A'/2A M/z'^, a' = da/dz ~ —2M/z'^ -C —a^, Ri^ujim = 

\{Rlt — — a') ~ M/z^, so p ~ M/z^, and then choosing Qr = 0 gives ur = 
-A-^/^Qz ~ -^z ~ ^'r = duR/dz QJ/z"^ > \aujR\, so R^^^t = ^r + 

aujR K, QJ/z‘^, and q 3 MJ/z^. Then Eqs. 


and (HZI) give 


ujq k, —?,J/z^ 


and 


cu± ^ ±^M/z^ + {3/2)J/z^ zs 


t'UJ't'UJ * 


(177) 


(178) 


Inserting Eq. ( [177|) into Eq. (|167|) then gives hlo ~ which agrees with the 

general far-held Eq. (|102|) with r = z. 

One can iterate the approximations of Eqs. (|175|) and (|176|) by rewriting the 


cubic Eq. (|169|) as 


Wo — 


-g + wg 

p 


(179) 
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and as 


LiJ± = it 



(180) 


Then one starts with setting the terms involving ujq or a;± eqnal to zero on the right 
hand sides, evaluates the right hand sides to get the first approximations for the left 
hand sides, enters these approximations back into the right hand sides to get the 
second approximations, and iterates to get the desired accuracy. This procedure is 
often a better way to proceed when -C (particularly when one has functional 
expressions for p and q in terms of some coordinate like z along the axis) than to 
use the explicit solutions (|173|) and (|174|) of the cubic Eq. (|169D . 


The discriminant of the cubic Eq. ( p.69|) , which is proportional to 


108g2 - 16p3 = - IQ{R 


2\3 


Hzutzu 


a) 


(181) 


is negative when the magnitude of the argument of the inverse sine in the exact 
solutions (|173|) and (|174|) , (27g^/4p^)^/^, is less than unity, leading to three real 
solutions for u. When one gets sufficiently deep into a strong rotating gravitational 
field that 27q^ exceeds 4p^, the roots ujq and a;_ merge at —^Jp/3 and then go off into 
the complex plane, leaving no real cuq for a SCAM but only the single remaining real 
root for the SCALE that has a global minimum of the acceleration, with normalized 
rotation 




4p3 


' 1 - 


4p3 


o'q. 


(182) 


8 The SCAM and other SCALEs 
in the Kerr-Newman metric 

Now let us evaluate the properties discussed above of the SCAM (Stationary Con¬ 
gruence Accelerating Maximally) and other SCALEs (Stationary Congruences Ac¬ 
celerating Locally Extremely) in the Kerr-Newman metric 

ds^ = —-[dt — a{l — c^)dp]‘^-\ - :^[{r^ + a^)dp — adt]‘^ + ^dr‘^-\ - - —(183) 

A 1 — 

where, to keep everything algebraic, I have used c = cos 6 instead of 6 in what are 
otherwise Boyer-Lindquist coordinates. Here, as usual (see, e.g., iSiiil) 

a = J/M (184) 

is such a standard Kerr parameter that I shall continue to use it in expressions 
directly involving the metric, even though elsewhere I use it for the acceleration, 

A = - 2Mr + (185) 
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which is D/{I — c^) in terms of D = gQ^ — 5'oo5'ii dehned by (^) and 


2 _ 2 I 2 2 

p = r + a c. 


The Kerr-Newman metric thus has 


, , 2Mr-g2 

A = —gtt — 1-5-— 1 — 2f/, 

p2 

a(2Mr — — c^) 2 \ 

B = -gt^ = ^^ = 2aU{l - c^), 


above, and 
(186) 

(187) 

(188) 


2 , „ ,, a^(mr - Q^)(l - cy 


C = -9,^ = -(r" + a^)(l-c-=) 

= —(r^ + a^)(l — c^) — 2?7a^(l — c^)^, 




(189) 


where 

_ 2Mr -Q^ _ Mr- Q'^/2 

2p2 J.2 _|_ ^ 2^2 

is a particular generalization of the Newtonian potential (with the sign reversed to 
make it positive for r > Q'^/ (2M)). In fact, another way to write the Kerr-Newman 
metric above is 



ds^ = -dt^ + (r^ a^)(1 - c^)dp^ + 

-I- 2f/{[(it — a(l — -I- 


2\j 2 , 7 2 


2 I 2 1 2 

(r^ -I- a^c^Ydr"^ 


-dc 


(r2 _|_ a2^^j,2 _|_ q2 _ 2Mr -|- Q"^) 

where the hrst line is simply flat spacetime in spheroidal coordinates P 

r = ^ — a?y + Aa^z "^], 

c = cos 9 = z/r. 


}, (191) 


(192) 

(193) 


One can in principle obtain the angular velocity O of the SCAM in the Kerr- 
Newman metric as a function of the coordinates r and c by evaluating Eq. (P^ and 
setting it equal to zero, which gives a quartic equation for O. However, when this 
equation is rationalized and expanded out in powers of a, M, Q, r, c, and O, one 
gets literally hundreds of terms. Thus it seems likely that the explicit solution would 
take more space to print than the entire rest of this paper, so I have not bothered 
to do that. 

In this way a SCAM, though it is simple to specify implicitly, is not nearly so 
simple to specify explicitly (e.g., by ff(r, c)) as a Zero Angular Momentum Observer 
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(ZAMO 0, 0) or even as an Extremely Accelerated Observer (EAO |^, ^), which 
is a Stationary Observer whose angular velocity extremizes the cylindrical radial 
component of the acceleration, and which can be given by a seven-line explicit 
expression for the Kerr metric (Q = 0). 

However, one can specify explicitly the SCAM in the three limiting cases de¬ 
scribed above for a general stationary axisymmetric metric. 


8.1 In the slow-rotation limit 


First, consider the case in which the Kerr parameter a is much smaller than M. This 


is the slow-rotation limit, and evaluating Eq. (|101|) to first order in a (but including 
the lowest-order term in 1/r which is cubic in the Kerr parameter a = J/M, the 
first term that depends on the angular variable c = cos 6) gives the angular velocity 
of the SCAM as 


On — 


-a{Mr^ — Q^r + 3a^c^) /a^M^ afQ 


^3(^2 _ 2,Mr + 2g2) 


A- 


aM r 3M2 

- T- 1 + 


This requires that 


Mr 


r > 


+ 


9M^ - 5Q^ + 3a^c 


/ 

- + ^( 


a^Q\ 


a^M 


r3 ’ ]\/[f3 /. 


;(3M+^9M2-8g2) 


.(194) 


(195) 


in order that the denominator of the first term on the right hand side of Eq. 
not have changed sign as one comes in from infinity, but the denominator may be 
arbitrarily small, so long as a is sufficiently small (e.g., not only small compared 
with M, but also much smaller than the square root of the denominator divided by 
Mr2). 

This has the consequence that, in the slow-rotation limit, the SCAM for Kerr- 
Newman is defined for values of the radial variable r obeying the inequality (|195|) , 
and its angular velocity is given by Eq. (|194|) to good accuracy so long as the 
magnitude of this expression for Hq is much smaller than 1/a. 

By comparison, in the Kerr-Newman metric the ZAMO has angular velocity 
given explicitly by 


Viz — 


a{2Mr - 


a{2Mr - 


(r^ -|- q2^2 _ 


0{a^). 


(196) 


Obviously, in the far-field limit, -|- g^, Eqs. 


and (^) agree 

with Eqs. ( p.02|) and (|103|) . However, outside the far-field limit but still within the 
slow-motion limit, one does not have the simple relation Viz = —2f2o that one has 
in the far-field limit. 
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8.2 At the orbitals (on the equatorial plane) 


Second, consider the case of the orbitals (locations of pairs of stationary geodesic 
observers, corotating and connterrotating tinielike circular Keplerian orbits). These 
all occur on the equatorial plane (c = 0) of the Kerr-Newman metric. There the 
angular velocities of the speed of light, given in general by Eq. (|35|) , are 


— 


±r^-\/r^ — 2Mr + + Q‘^ + a(2Mr — Q^) 


(197) 


j,2(^2 _|_ gj2^ _|_ (P{2Mr — Q^) 

and the angular velocities of the SCALEs that are the Keplerian orbiting stationary 
geodesics, given in general by Eq. (|109|) , are 


±y/Mr - Q'^ 
± a^/Mr — 


(198) 


On the equatorial plane the other SCALE, the SCAM which gives a local max¬ 
imum of the acceleration, reduces to what Semerak |]3^ , ^ calls an Extremely 
Accelerated Observer (EAO), with Eq. (| 111 |) giving its angular velocity as 


r2(r2-3Mr + 2 Q 2 )_ 2 a 2 (Mr- + 

2a[r2(3Mr - 2Q‘^) + a?{Mr - Q^)] 

_ -2a{Mr - Q^) _ 

r^{r^-3Mr + 2Q^)-2a^{Mr-Q‘^)+r^^{r^-3Mr + 2Q^y-Aa‘^{Mr-Q^) 
-a{Mr - Q^) a^{Mr - Q‘^){2r‘^ - 3Mr + 2Q‘^) 4 

r2(r2-3Mr + Q2)^ + ^ 2(^2 _ + 2Q^y ^ 

a, 3 M 2 -g 2 9M3-5MQ2 27M^-21MQ‘^ + 2Q^ + 2a^M^ 

-- M+-^ ++-+0 - , 

(199) 


a straightforward extension to Q 7 ^ 0 of the result in the Kerr metric 

The velocities of the circular Keplerian orbits in the SCAM frame (both of equal 
magnitudes but of opposite signs, as Semerak found was the case in the Kerr equa¬ 
torial plane and later found in general [Q) have magnitude given by Eq. (|117|) 
specialized to the Kerr-Newman metric: 


Vk = 


r(r - M) - ^(r2 - 3Mr + Q 2)2 _ 4^a?{^Mr - Q^) 


2^{Mr - Q2)(r2 - 2Mr + a2 + Q 2 ) 

2^{Mr - g 2 )(r 2 _ 2Mr + a2 + g 2 ) 
r(r - M) + ^(r2 - 3Mr + g2)2 - 4a2(Mr - g2) 

a?r{r — M) 


Mr - g2 


■[1 + 


r2-2Mr + g2' 2{r‘^ - 2Mr + Q‘^){r‘^ - 3Mr + 2Q‘^f 


+ 0{a% 


( 200 ) 
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Similarly, Eq. 


specialized to the Kerr-Newman metric gives the magnitude of 


the acceleration of the equatorial SCAM as 

VMr - Q'^ 


Oq — 


-Vk 


r(r - M) - - 3Mr + - Aa^{Mr - Q^) 

2r‘^yjr‘^ — 2Mr + + Q‘^ 

2(Mr - Q^)y/r^ - 2Mr + a'^ + Q'^ 


r2[r(r - M) + - 3Mr + Q 2)2 _ 4a^(^Mr - Q^)] 

Mr — ^ a?‘r{r — M) 


^[1 + 


- 2Mr + ^ 2{r^ - 2Mr + Q^){r^ - 3Mr + 2Q^y 


+ 0(a")] 


( 201 ) 


Then one can use Eq. (|127|) to get the acceleration of a stationary observer at any 
speed V relative to that of the SCAM on the equatorial plane. 

Many of these formulas for the properties of the SCAM on the equatorial plane 
of the Kerr-Newman metric are simpler in the case of extreme Kerr-Newman, = 
— a^, so that A = {r — MY and the event horizon is at r = M. Then, taking a to 
be positive, one can calculate that the SCAM exists for values of r down to 2M-|-2a, 
where a dehned by Eq. ( |119|) passes through zero and goes negative, making flo and 
various other quantities complex for smaller r. Hence I shall hrst give the value at 
general r when = M^ — a^, and then, after the hrst arrow in each equation, the 
limiting value for each quantity when one sets r = 2M -|- 2a. Finally, after a second 
arrow, I shall give the limiting value in extreme Kerr a = M (so Q = 0) at r = 4M, 
the radial inner boundary on the equatorial plane of the region where the SCAM 
exists as a real congruence: 

r‘^[r — M) + a{2Mr — M‘^ + a?) 


f2_i_i — 


M -|- 3a 

r2(r2+a2)+a2(2Mr-M2 + a2) ^ 4M^+9aM+7^ 


= 


-r2(r-M) + a(2Mr-M2 + a2) 
r2(r2 + a2) + a2(2Mr-M2 + a2) 

VMr -M^ + a? 


-1 


^+K — 


AM + 3a 
1 


r2 -1- a\/Mr — M^ -|- a^ 


Q.-K — 


-YMr - M2 + 


4M -F 5a 
-1 


9M’ 

1 


r2 — aylvP + a? AM + 3a 


7M’ 


1 

5M’ 

(202) 

1 

7M’ 

(203) 


(204) 


(205) 


Hn — 


-2a(Mr - M^ + a^) 


r2 (r - M) (r - 2M) - 2a2 (r2 + Mr - M2 + a2) + r2 (r - M) ^(r - 2M)2 -4a2) 

-11 . ^ 

(206) 


4M -|- 3a 


7M’ 
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1 




a{2Mr - + a^) 


a(3M + a) 


Vk 


{r‘^ — ar + aM+a‘^){r‘^ + ar — aM+a‘^) 
r — \J {r — 2M — 2a){r — 2M+2a) 


(4M+3a)(4M2 + 9aM+7a2) 
2VMr -M^ + a? 


2VMr + r+^{r-2M-2a){r-2M+2a) 


ao — 


r — 2M — 2a) (r — 2M + 2a) 


35M’ 

(207) 

^ 1 , 

(208) 
(209) 


2r2 4(M + a) 8M’ 

One can thus see that at the inner boundary of the SCAM, at r = 2(M + 
a) on the equatorial plane of extreme Kerr-Newman (Q^ = — a^), both the 

counterrotating circular Keplerian orbit and the SCAM have angular velocities, ^-k 
and Oo respectively, that approach the angular velocity 0_i of the counterrotating 
speed of light. Since at this inner boundary = 1, Eq. (|127|) says that the 
acceleration there is independent of the (signed) speed n of a stationary observer 
relative to the SCAM if |n| < 1. However, one must take care, since this formula 
is then degenerate at |n| = 1, and since Eq. (|116|) with Hi = Hq = gives 
n(H2,Hi) = 1 for any angular velocity H2 > Hi. For H > Hq = H_i, one should 
instead return to Eq. ([31|), which gives at r = 2(M + a) in extreme Kerr-Newman 
the acceleration 

(M + 2a)[l-(4M + 5a)H] 3 /1-9MH\ 


a = 


4(M + a) [(M + 3a) - (4M2 + 9aM + 7a2)H] 32M V1 




( 210 ) 

with the expression after the arrow being that at r = AM when a = M (and hence 
Q = 0). Note that Eq. (|210|) is one where the a on the left hand side is the 
acceleration, whereas all the a’s on the right hand side denote the Kerr parameter 
a = J/M^ as I have warned. 

8.3 On the axis of rotation 

The third limiting case where one can give the SCAM explicitly for the Kerr-Newman 
metric without solving a very messy quartic equation is on one of the axes of sym¬ 
metry, say, for concreteness, the one at 0 = 0 (c = cos 6* = 1). There the rotation ujq 


of the SCAM is given by the solution (|173|) of the cubic Eq. (|169|) , uj^ — pu — q = 0, 


where the coefficients p and q, given by Eqs. ( |170| ) and (|171| ) in the general case, 
take on, in the Kerr-Newman metric, the values 

Mr^ — Q'^r'^ — Sa'^Mr + — Q'^r — a^M)’^ 


p = 


and 


(r2 -|- a2)3 (r2 -|- — 2Mr -1- a2 -|- 

a{Mr^ - Q^r - a^M){3Mr‘^ - 2Q^r - a^M) 


q = 


(r2 -1- a2)9/2(r2 — 2Mr -|- a2 -1- 


( 211 ) 


( 212 ) 
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Once one has the normalized rotation which is independent of the choice 
of the Killing vector field k (though near the axis 1 is uniquely determined, up to 
sign, by its property of having a magnitude which goes as the proper distance from 
the axis), one can calculate, for k = d/dt, the angular velocity Oq of the SCAM 
arbitrarily near the axis by solving Eq. (|167|) . On the axis of the Kerr-Newman 
metric, this gives 


Oq — + \/~Au]q 

q{2Mr - 


^ I 


+ a?) 




2Mr + a2 + Q2 
r2 + a2 


- sm - sm 
3 ^3 


-1 


27 




In the case of the extreme Kerr metric, a = M and Q = 0, if we let x = r/M 
and r/ = (r^ + /M'^, then the cubic Eq. ( [I69|) becomes 

— {x^ — — bx — l)y — (3a;^ + 3x^ — a; — 1) = 0. (214) 

The discriminant of this cubic is proportional to 31 + 114a;+115x^ + 56a;^+12a;‘^—4a;^, 
which is negative for 


x = — 6.15862999016071260705, 

M 


(215) 


the region where the cubic has three real roots (one for the SCAM and one for each of 
the other two SCALEs that are local minima of the acceleration). Therefore, near the 
axis of an extreme Kerr black hole, the SCAM exists only outside r ^ 6.15862999M, 
whereas on the equatorial plane it exists outside r = AM. 

Although it is rather messy to do for general a and Q, for extreme Kerr one can 
also readily compare expansions of hlo iii inverse powers of r in the equatorial plane 
and on the axis. In the equatorial plane for a = M and (5 = 0 one gets 

3M 9M2 29M3 

fio =- ^[1 +-+ ^ + —^ + 0 (^)], 


whereas on the axis one gets 

M2 


^ ... , 3M 12M2 53M3 ^,M\, 

= + + +0(^)]. 


(216) 


(217) 


The difference in the two expressions, starting at the second-order correction, per¬ 
sists even when one changes the radial variable from r to p = \/r‘^ + a? cos2 9 (which 
is the same as r on the axis), since then the series on the equatorial plane becomes 

3M 21M2 35M3 


fin — 


M2 

p3 


[1 + 


P 


+ 


2p2 


+ 




+ 0(^)], 


(218) 


whereas on the axis it has the same form as Eq. (pi7]) but with r replaced by p. 
The first three terms of Eqs. (|216|) and (|217|) can be readily be seen to agree with 


Eq. (|194|) when a = M, (5 = 0, and either c = 0 (0 = 7r/2) or c = 1 (6* = 0). 


46 
































9 Stationary Worldlines Accelerating 
Radially Maximally (SWARM) 

Because the explicit expression for the angular velocity hlo of the SCAM (a root 
of a messy quartic equation) appears to be generally rather intractable except in 
the special cases discussed above (slow rotation, at an orbital, and on an axis), it 
might be useful to define other congruences that have simpler explicit expressions 
and which have some of the properties of the SCAM. For example, one might define 
Stationary Worldlines Accelerating Radially Maximally (SWARM) as those that 
have n chosen to maximize the radial component of the acceleration. 

Then the question arises as to how to define the radial component (or direction, 
since the radial component is needed only up to a positive constant of proportionality 
at each location in order to be able to define its maximum as a function of the angular 
velocity). Semerak [^, defined a maximum of the cylindrical radial component of 
the acceleration in the Kerr metric as an Extremely Accelerated Observer (EAO), 
but I prefer to dehne something different here and leave that name for what he 
dehned there. (And to avoid confusion I propose that the EAO retain its original 
definition in [^, rather than being redefined to be the SCAM that I invented, 
despite Semerak’s proposal to do that in |l3^ .) 

Alternatively, if one dehned the radial direction to be that of the acceleration of 
the SCAM, then of course the SWARM would simply be the SCAM, but then Ending 
the radial direction would involve solving a quartic equation, and there would be no 
advantage to dehning a SWARM. 

Therefore, I propose that a SWARM be dehned so that in the Kerr-Newman 
metric the radial direction is that of Vr, the gradient of the Boyer-Lindquist radial 
coordinate r. This direction has several remarkable properties in Kerr-Newman, 
connected with the existence of Carter’s ‘fourth constant of motion’ 


37, 


For example, there is the fact that for a timelike geodesic, out of the set of four 
parameters that govern the orbit in the (2,3) or (r, c = cos 6*) coordinates (e.g., the 
value of r at c = 0, the value of dr/dc there, the value of the conserved energy 
per rest mass —Uq = —Ut, and the value of the conserved angular momentum per 
rest mass ui = u^), a two-parameter subset leads to orbits with constant r (e.g., r 
and Ml, setting dr/dc = 0 at c = 0 and choosing —uq as a function of ui so that 
d?r/dc^ = 0 there). For a generic stationary axisymmetric metric, one could dehne 
a radial coordinate so that a one-parameter set of geodesics have constant r (e.g., r 
along some one-dimensional line in the (2,3) plane analogous to the c = 0 line, by 
dehning r so that it is constant along the orbit that one gets with a specihc choice of 
—Mo and Ml), but if one tried varying a second parameter, e.g. mi, then there would 
be no choice of —Mq for diherent values of Mi that would give other orbits along the 
same constant r line. 
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In the Kerr-Newman metric, the geodesic orbits with constant r typically oscil¬ 
late in c (unless they are circular Keplerian orbits in the equatorial plane, which 
stay at hxed c = 0 as well as constant r and hence are stationary geodesics), go¬ 
ing between a positive maximum where dc/dr = 0 as well as dr/dr = 0 (so that 
the velocity is momentarily zero in the (2,3) plane at that turning point, though of 
course the velocity is not zero in the (0,1) plane), and a negative minimum of equal 
magnitude (because of the symmetry with respect to reflections about the equatorial 
plane) where also the velocity in the (2,3) plane is momentarily zero, dc/dr = 0 as 
well as dr/dr = 0. 

However, in a generic stationary axisymmetric metric, an orbit that starts with 
zero velocity somewhere in the (2,3) plane will not generically have zero velocity 
elsewhere on its orbit in that plane. Nevertheless, requiring zero velocity in the 
(2,3) plane at one point there puts only one restriction on the conserved quantities 
—Uq and Ml, so by choosing wi, say, appropriately (which then leaves —Uq determined 
as a function of ui so that the velocity in the (2,3) plane is zero there), one has the 
right number of free parameters to be able to get zero velocity in the (2,3) plane at 
some other location. If this requirement of a second turning point of zero velocity 
in the (2,3) plane uniquely hxes ui (and hence also —mq) at the original turning 
point, then the direction in which the orbit starts out (determined by d?x°'/dT'^ at 
the turning point where dx°‘/dT = 0) could be dehned as the angular direction, 
the direction of constant radius r (after dehning the radius r suitably). Then the 
orthogonal direction could be dehned as the radial direction. 

This is admittedly a nonlocal dehnition of the radial direction, and it might not 
always give a unique answer, but it is a procedure that would give the Vr direction 
in the Kerr-Newman metric and presumably would give a unique direction for small 
perturbations of that metric. So let me dehne the radial direction by this procedure 
when it works, and then dehne a prime as denoting a partial derivative in that 
direction. 

Then Eq. (|1 1 1|) gives the angular velocity of what is now the SWARM, which 
maximizes the radial component (rather than the entire magnitude) of the acceler¬ 
ation. Eq. ( |109| ) gives the angular velocities, not of circular Keplerian orbits, but of 
stationary worldlines with zero radial component of the acceleration. These will also 
be the angular velocities at that location of the nonstationary geodesic orbits that 
have a turning point in the (2,3) plane at that location and also have d'^r/dr'^ = 0 
there. By the same argument as before, these have equal but opposite velocities 
in the frame of the SWARM, of magnitude given by Eq. (|117|) . Furthermore, Eq. 
( 125 ) gives, not the total acceleration, but the radial component of the acceleration 
of the SWARM at that location, and Eq. (|127|) gives the part of the acceleration in 
the radial direction for a stationary worldline at speed v in the SWARM frame. 

In the Kerr-Newman metric (|183|), one can write out all these expressions explic- 
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itly as functions of r and c, but they don’t easily fit into single lines, so I won’t bother 
doing that straightforward calculation here. One can see that in the Kerr-Newman 
metric, which has a two-parameter family of timelike geodesics with constant r as 
mentioned above, the SWARM four-velocity is the normalized average of the four- 
velocities of the two constant-r geodesics that have their turning points in the (2,3) 
plane at that location (i.e, have their maximum of |c| there). In other words, if one 
takes two equal-mass particles moving along constant-r geodesics that at a certain 
location in the (2,3) plane have their four-velocity entirely in the (0,1) plane, with 
one moving forward in ip and the other moving backward, then if one makes a to¬ 
tally inelastic collision between these particles, immediately after the collision the 
resulting particle will have the four-velocity of the SWARM at that location. 

One can also say that on the equatorial plane in the Kerr-Newman metric, as at 
any orbital in a generic stationary axisymmetric metric, the total acceleration is in 
the radial direction, so the SWARM has the same four-velocity as the SCAM there. 
The four-velocity (and acceleration) also agrees on an axis, trivially, since there the 
four-velocity is independent of the angular velocity, but the angular velocity of the 
SWARM is not the same as that of the SCAM near the axis of the Kerr-Newman 
metric. In the limit of going onto the axis itself, Eq. ( pi3|) gives the angular velocity 
of the SCAM, and the angular velocity of the SWARM there is 


ffn — 


—a(Mr^ — Q'^r — a?M) 


(r^ -|- a?){r^ — 3Mr -|- 2Q‘^r -|- a^r -|- a?M) 

9M2 - 5g2 _ 2a2 


aM r 3M2 
1 + 


Mr 




a^Q\ 


a^M 


r3 ’ 


(219) 


By comparing with the series expansion of Eq. ( [194| ) at the axis (c^ = 1), one 
sees that, to lowest order in 1/r, the angular velocity of the SWARM on the axis is 
less negative than that of the SCAM by 5a^M/r^, so if a/M is small, the SWARM 
is a very good approximation for the SCAM on the axis, and, presumably, at all 
other angles or values of c. This agreement is a consequence of the fact, that for 
small a/M at least, the direction of the acceleration of the SCAM is nearly radial in 
the Kerr-Newman metric, so maximizing the radial component of the acceleration 
for the SWARM (Stationary Congruence Accelerating Radially Maximally) gives 
very nearly the same angular velocity as finding the local maximum of the magni¬ 
tude of the total acceleration for the SCAM (Stationary Congruence Accelerating 
Maximally). 


10 Application to Maximally Rotating Stars 

One application of the fact that the SCAM is usually counterrotating relative to a 
rotating source (e.g., a black hole or star) is an explanation of the fact that corotating 
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Keplerian orbits in an equatorial plane of the source may have periods that are longer 
than one would get from a naive application of Kepler’s third law. In certain cases 
Kepler’s third law gives a better approximation to the orbital frequency hi — hlo 
relative to the SCAM, rather than the orbital frequency hi relative to a nonrotating 
observer with Qnr = 0. Since hlo is typically negative (when the coordinate system 
is chosen so that the source is rotating positively), the orbital frequency relative to 
a nonrotating observer will be slower than that relative to the SCAM, so its period 
will be greater. 

Physically, this effect may be explained by the same mechanism used above 
to explain the counterrotation of the SCAM: A corotating orbiting observer will 
see the part of the source nearest her as partially moving with her and hence as 
having a lower energy density in her frame than the part of the source farthest 
away from her, which is moving in the opposite direction. Hence she will see the 
energy distribution shifted slightly away from her, where it will have a weaker net 
gravitational attraction on her than a source of the same energy density distribution 
(as seen in a nonrotating frame) that is not rotating. Therefore, she will orbit more 
slowly by this relativistic effect that can be ascribed to the angular momentum of 
the source. 

Realistic situations are further complicated by the fact that the source (e.g., a 
star) will not generally have the same energy density distribution in a nonrotating 
frame when it is rotating as when it is nonrotating. One effect is that if a source 
is spun up, it will gain rotational energy and hence total mass-energy. However, 
this effect can easily be compensated for by removing from the source an amount 
of energy equal to that given it in spinning it up, say by removing an appropriate 
number of baryons. In any case, if one is using Kepler’s third law, the mass in that 
formula should be the total mass-energy in the source, so a change in the total mass 
is already taken into account. 

However, another effect that is not taken into account by Kepler’s third law is 
that the shape of the source generally changes when it is spun up. For example, 
if the source is a self-gravitating fluid, such as a star, the centrifugal forces of the 
rotation will generally cause the star to become oblate. Then as the energy becomes 
more concentrated upon the equatorial plane, where it is on average nearer to the 
observer orbiting in that plane, it will exert a greater gravitational attraction upon 
the observer, causing her to orbit faster than she would have in the absence of the 
oblateness. Equivalently, the greater gravitational attraction at a hxed radius in the 
equatorial plane is caused by the quadrupole moment of the source. 

This oblateness or quadrupole effect on the orbital frequency is opposite in sign 
to the relativistic effect of the source angular momentum discussed above. At low 
source angular velocities, the relativistic effect is linear in the source angular ve¬ 
locity, whereas the oblateness effect is quadratic. However, the quadrupole effect 
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persists even in the Newtonian limit, so it can be larger than the relativistic angular 
momentum effect. 

As an example in which the corotating orbital period in the equatorial plane 
can be calculated exactly and compared with Kepler’s third law, consider the Kerr 
metric with zero charge. Then from Eq. (|198|) , and temporally restoring Newton’s 
constant G and the speed of light c that have been set equal to unity, one can readily 
get the period as 


P = 


27r 


n 


= 2ti\ 


+K 


I 2'nJ 


( 220 ) 


The second term on the right hand side can be identihed with the linear increase in 
the period with the angular momentum J, and the square of the speed of light in 
the denominator shows that it is a relativistic effect. (It is interesting that Newton’s 
constant does not appear in this term. If one takes the quantum-mechanical phase 
of a system with energy E = Mc^ and angular momentum J to be 
then the time period for this phase to rotate around once is precisely this second 
term. However, Newton’s constant does appear in the ratio of the second term to 
the hrst term, so the increase in the period with the angular momentum really does 
involve both special relativity and gravity and is thus a general relativistic effect 
rather than purely a special relativistic effect.) 

It is tempting to identify the hrst of the two terms on the right hand side of 
Eq. ( | 22 U|) as being precisely Kepler’s third law for a circular orbit of radius r, but 
there is the question of whether r is the most natural measure of the the radius. In 
the nonrotating limit (Schwarzschild, J = Ma = 0), r is the circumferential radius, 
l/27r times the proper circumference of a closed circle in the equatorial plane with 
hxed r and t, which is a fairly simple geometric dehnition of a radius that makes 
Kepler’s third law exact for circular orbits in the Schwarzschild metric. If in this 
section one takes R to be the circumferential radius in the equatorial plane around 
a source of total mass M, then one can dehne 


GM 


( 221 ) 


as a measure of how closely Kepler’s third law holds, which would state that T = 1. 

Although T = 1 for the Schwarzschild metric, in the equatorial plane of the Kerr 
metric with J 7 ^ 0 the circumferential radius is 


so 


R = ^/^ 



J2 

M2c2 


2GJ^ 

d^Mr' 


( 222 ) 
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/ 2GJ^ 

^ McV3 
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= 1 - 


2J G 


+ 


(? V 2M'^c^r‘^ 


+ 0 (r-"). 


Alternatively, one can solve Eq. (|222|) for r as a function of R and insert this into 
Eq. (|22CI|) to get the orbital period (as seen from radial inhnity) as 


P = 


27r 


n 


= 2ti\ 




+ 


2'kJ 


37rj2 


+K 


GM Mc^ 2 M^c'^VGMR 


+ 0(i?-3/2). 


(224) 


The third term on the right hand side represents the effect of the quadrupole moment 
of the Kerr metric. 

Since in the Kerr metric a stationary observer must have r > GM/c^, and since 
the dimensionless Kerr rotation parameter, 


cJ 




M GM^' 


(225) 


is less than or equal to unity, the negative of the ratio of the second term to the 


third term in the last expression for T in Eq. ( P23|) , or in the last expression for P 
in Eq. (|2^ , is 

jlQGM^r _ 4 / c^r 

V 9J2 

Therefore, at least at large r where one can drop the 0(r~^) term in Eq. (|223|) 


> 1 . 


(226) 


or the 0{R~^^^) term in Eq. (P24|) , the second term of either of these equations, 
which represents the relativistic period-increasing effect that is linear in the angular 
momentum J, dominates over the third term, which represents the effect of the 
period-decreasing quadrupole moment that is quadratic in the angular momentum. 

However, if the dimensionless Kerr rotation parameter a is larger than about 
0.952518, then there are stable circular corotating orbits in Kerr (which exist for 


— QMr + 8aVMr — 3a^ > 1, 


(227) 


temporarily reverting to units in which G = c = 1) for which the effect of the 
quadrupole moment dominates so that T > 1. For example, for the extreme Kerr 
metric (a = 1 or a = M), T > 1 for r less than about 2.01186GM/c^, and at 
the smallest innermost stable circular corotating orbit at r = GM/c^, one has 
R = 2GM/c^, Q+k = c^/2GM, and hence T = 2. 

If we turn to models of maximally rotating stars (stars with the equatorial sur¬ 
face rotating at the Keplerian velocity and hence just marginally bound), typically 
the quadrupole moments are larger than they are for the Kerr metric with the same 
stellar mass and angular momentum, because stars are not so gravitationally con¬ 
centrated. Thus one gets a larger radius at which the relativistic period-increasing 
effect (linear in the angular momentum) balances the effect of the period-decreasing 
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quadrupole moment. For nonrelativistic maximally rotating stars this radius tends 
to be outside the surface of the star, so that one usually gets F > 1 from Eq. (|221|) 
when R is set to be the circumferential equatorial radius of the star in its rotating 
frame, and VLj^k is both the angular velocity of the corotating circular equatorial 
orbit at the surface of the star, and the angular velocity of the star itself. 

However, for certain relativistic maximally rotating stars, the relativistic period- 
increasing effect can exceed the period-decreasing effect of the quadrupole moment, 
making F < 1. For example, in the numerical models of rapidly rotating polytropes 
in general relativity by Cook, Shapiro, and Teukolsky |^|, the maximum uniform 
rotation models of the “supramassive” sequence given in their Table 2 have F < 1 
for values of the polytropic index n < 1.5, namely F = 0.958 for n = 0.5, F = 0.988 
for n = 1.0, and F = 0.999 for n = 1.5. Higher values of the polytropic index, 
n > 2, generally seem to give F > 1, namely F = 1.005 for n = 2.0 and F = 1.009 
for n = 2.5. The data from 1^^ for n = 2.9 naively gives F = 0.999 < 1, but since 
the data are given only to three places, and since this large value of n gives a highly 
nonrelativistic model {2GM/Rc^ = 0.0109, as opposed to the highly relativistic value 
2GM/Rc^ = 0.578 for n = 0.5, for example), I would be sceptical that actually F < 1 
for this large value of n. 

Similarly, one might doubt that F < 1 for n = 1.5, where the three-place data 
give a result also just slightly below unity, but presumably there is a value of the 
polytropic index n fairly near 1.5 such that F = R?Vl\j^/[GM) < 1 for maximally 
uniformly rotating supramassive polytropic models with smaller n but such that 
F > 1 for similar models with larger n. The small values of n give relativistic stellar 
models for which the counterrotation of the SCAM is sufficient to make it so that 
the corotating Keplerian orbital velocity (which matches the stellar rotation rate at 
the equatorial surface in these maximally rotating models) relative to a nonrotating 
observer at infinity is slowed down by this relativistic effect more than it is sped up 
by the quadrupole moment. 


11 Application to the gravitational field 
of the Sun and Solar System 


Another application of the stationary congruences defined above and of the related 
deviations from Kepler’s third law is to the gravitational field of the Sun. Since 
the Sun is nearly spherical and is not very relativistic, the metric for its external 
gravitational field may be given, to an accuracy of about one part in 10^®, and using 


units in which (7 = c = 1, as ||40 


ds^ = 


[1 


2Mr: 


© _l_ 2Q0 


P2{cos9)]dt^ + [1 


2Mq 2Qq 

-^- —P2{cos 9)\ dr 
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(228) 


_l_ [X _ ?^®P2(cos6')]r^[(i6'^ + sin^ 9{dip — 

y^O y»0 

Here P2{cos9) = {3cos^ 9 — l)/2 is the standard second-order Legendre polynomial, 
and Qq is the qnadrnpole moment of the Snn (not to be confnsed with the previons 
nse of Q to denote the charge of the Kerr-Newman black hole; in this section the 
charge will always taken to be zero and Q will always denote a qnadrnpole moment, 
with its sign chosen so that an oblate spheroid, snch as a rotating body, has positive 

Q), 

Qq = J 2 MqRq = — J P 2 {cos 9)dr sin 9d9dip = J p( 2 ^^ + 2 ^^ ~ z^)dxdydz. 

(229) 

Here Rq is the radins of the Snn, and J 2 is its dimensionless qnadrnpole moment 
parameter. 

In metric length nnits (meters), the 1996 Review of Particle Physics gives 
the mass and radins of the Snn as 

GMq/c^ = 1476.625 04 m = 0.9137 x 10^® ip, (230) 

= 6.96 X 10® m = 4.31 X 10^^ £p, (231) 

where ip = \JhG/c^ = (1.616 05 ± 0.000 10) x 10“®® m is the Planck length, nsing 
data from the same sonrce ||41|| . 

A recent helioseismic determination of the solar gravitational qnadrnpole mo¬ 
ment |42|, which is consistent with a less-precise direct measnrement of the solar 
oblateness [^], gives the dimensionless solar qnadrnpole parameter as 

J 2 = (2.18 ±0.06) X 10“^. 

Therefore, the solar qnadrnpole moment is 

Qq = J2MqR% = (2.10 ± 0.06) X 10^^ kg m^ 

or in length nnits (cnbic meters) it is 


(232) 


(233) 


GQq/P = GJ 2 MqR^q/P = (1.60 ± 0.04) x 10^^ m® = 160 000 ± 4 000 km 


= (53.8 ± 0.5 km)® = (3.69 ± 0.11) x 10^^^ i],. 

From this one can define an effective qnadrnpole radius of the Sun as 


(234) 


tqq = ^J2Qq/Mq = 460 ± 6 km 


= (3.07 ± 0.04) X 10“® AU = (2.84 ± 0.04) x 10^° ip, 


(235) 


the radius of a solar-mass ring of the same qnadrnpole moment as the Sun. Here 
1 AU = 149 597 870 660 ± 20 m = 0.9257 x 10"^® ip is the astronomical unit O . 
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The same helioseismic measurements also give the angular momentum of 
the Sun as 


Jo = (1.900 ± 0.015) X 10^^ kg m^ s“^ = (1.801 ± 0.014) x 10^^ h, 

which when converted to length and area units gives 

GJq/c^ = (4.705 ± 0.037) x 10^ m^ = (686 ± 3 m)^ 

= 47.05 ± 0.37 hectares = 116 ± 1 acres. 


(236) 


(237) 


From this and the solar mass one can readily calculate that the Kerr parameter a of 
Eq. (|184|) and the dimensionless Kerr rotation parameter a of Eq. (|225|) have the 


values 


a© = 


Jq 


= 318.65 ± 2.52 m = (1.972 ± 0.016) x 10^^ ip, 


MqC 


_ QqC^ _ cJq 
a© = 


GMp 


0 


GMl 


= 0.2158 ±0.0017. 


(238) 


(239) 


The fact that a© < 1 means that if the Sun were able to undergo gravitational 
collapse to become a black hole (which it is not, since it is too light, except for 
some extremely tiny tunneling probability or possibly some artificial compressional 
procedure), it would not need to lose any angular momentum to do so. 

By comparing the metric (|228|) with the uncharged Kerr metric (|183|) , one may 
deduce that the quadrupole moment of the Kerr metric is Q = J^/Mc^ = Ma^, 
so Mc^Q/J'^ = 1. (Remember that in this Section, Q denotes a quadrupole moment 
and not the charge, which we are here setting to zero.) For the Sun we find that 


Qe 


0 


J|/M©, 


= (1.04 ±0.05) X 10^ 


(240) 


so the quadrupole moment of the Sun is about a million times larger than that of a 
Kerr metric with the same mass and angular momentum. This is reasonable, since 
the dimensionless quadrupole moment J 2 is two-thirds of the solar oblateness WM, 
which one expects to be of the order of the square of the angular velocity of the Sun 
divided by the Kepler orbital velocity at the surface of the Sun, which is 
QqRq/GMq. Therefore, one expects (to order of magnitude) Qq = J 2 MqRq ~ 
RqQq/G. Then since J© ~ MqR^Qq, one gets M©c^(5©/J© ~ RqC^/GMq = 
0.471 X 10 ®, within about a factor of two of the correct answer. 

From the angular momentum J© of the Sun, one can calculate that the Stationary 
Congruence Accelerating Maximally (SCAM) in the equatorial plane rotates around 
the Sun with an angular velocity (as seen by a nonrotating observer at infinity) 

^0 ~ = -( 4.213 ± 0 . 033 ) X 10-20 (—) s"^ 

\ r J 
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_27r_ 

476 000 ± 4 000 yr V r y ’ 

and a linear velocity 


T / ATI\^ 

no ^ ^ = (6.303 ±0.050) X 10“M-) ms"^ 

V r / 

= -(2.932 ± 0.023) X 10"M—) ms"^ 

/ R \ ^ 

= -(9.253 ±0.073) km/yr. 


(241) 


(242) 


This rotation rate is extremely slow, thongli not qnite glacially slow, bnt in the 
4.6 ± 0.1 billion year age of the Solar System, the member of the SCAM at the 
snrface of the Snn would have made almost ten thousand backward revolutions 
relative to the distant stars, assuming that the far-held Eq. (|102|) applies all the 
way down to the surface of the Sun. 

Since the SCAM maximizes the magnitude of the acceleration for all stationary 
observers at a given location, it has a larger acceleration than that of a nonrotating 
observer. For the weak held of the Sun, the diherence is very tiny: 


Aa 


ao aNR ^ ^ ^ 


(2.656 ± 0.042) x 10"^® 

■AU\® 

- m s 

r / 

(4.775 ± 0.075) x 10"®° 

, r / 


(243) 


Over the lifetime of the Solar System, this tiny acceleration at r = 1 AU would, 
starting from rest, give a spatial motion lAat^ of about 2800 km. 

If, despite its rotation, the Sun were perfectly spherical and hence had no 
quadrupole moment, then Kepler’s third law would apply to high accuracy to the 
orbital angular velocity relative to the SCAM. Since the SCAM is counterrotating, 
the corotating Keplerian orbit has a lower angular frequency, and hence a longer 
period, relative to a static observer than relative to the SCAM. However, for the 
realistic Sun, the quadrupole moment from the oblateness, also caused by the ro¬ 
tation, increases the equatorial gravitational attraction, and hence also the orbital 
angular frequency there. This effect thus decreases the orbital period as measured 
by a nonrotating observer at inhnity. These two changes in the period may be 
calculated as follows: 
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In the metric Eq. 


for the gravitational held of the Snn, in the equatorial 


plane {6 = 7r/2 or P 2 {cos 6 ) = —1/2) the circumferential radius is 

Qq 




Qjp 


r + 


2 r 2 ' 


(244) 


The corotating Keplerian orbital velocity is then 

^ Mq ( Jq 


^+K — 
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3Qq 


+ OiR 


-3\ 


sJMqR^ 4:MqR^ 

and so the orbital period (as seen by a nonrotating observer at inhnity) is 
271 


(245) 


P = 
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+K 


= PK + APj^+APj, + OiR-^/^) 


= 271 \ 


R^ 
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277 J, 


© 




(246) 


GMq MqC^ 2^GMlR 
where Newton’s gravitational constant G and the speed of light c have been restored 


in Eq. (|246|) . 
Here 


Pk — 2,71 \ 


I R^ 

GMq 


= 31558 196.0 


R 


(m) 


(247) 


is Kepler’s third law for the period of a test body in the held of the Sun (ignoring 
the gravitational ehects of the planets for now), 


^Pje = 


271 Jq 
MqC^ 


= 277aQ/c = (2 002 ± 16 m)/c 


= (6.678 ± 0.053) x 10"® s = 6.678 ± 0.053 /is 


(248) 


is the increase in the period due to the linear ehect of the Sun’s angular momentum, 
and 

=-^ -( 1.117 ± 0.032) x 10-“ (A 

is the decrease in the period due to the quadrupole moment of the Sun. 

One can see that these corrections to Kepler’s third law in the gravitational 
held of the Sun are very small and currently unmeasurable, but it is amusing to 
calculate them as an academic exercise. It is also amusing to note that the ehect of 
the quadrupole moment (which is essentially quadratic in the angular momentum) 
dominates over that linear in the angular momentum for radii r < tkq, where 



tk 


0 


IQGJIMq 


280 ± 20 AU 


(250) 
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is the orbital radius at which Kepler’s third law would be exact, at about seven 
times the orbital radius (semimajor axis) of Pluto. 

One can use the order-of-magnitude estimates above for the angular momentum 
Jq MqR^VLq and the quadrupole moment Q© ~ to estimate 




Rq, 


(251) 


where Vr = Rq^q is th e linear rotation velocity of the equatorial surface of the Sun 
and Ve = ^J2 GMq/Rq = 2.060 x 10“^ c is the escape velocity from the surface of the 
Sun. The Sun is not rotating rigidly, so its angular velocity O© is not constant, but 
one can take, as a sort of averaged value for hi©, twice T©, the total kinetic energy 
in rotation of the Sun, divided by the Sun’s angular momentum J©. Since the total 
kinetic energy in rotation is ^ 


T© = (2.534 ± 0.072) x 10^® kg m' 




one gets an effective averaged angular velocity of the Sun as 

2T© ^—R _i 27r 


hi© — 


= (2.67 ±0.10) X 10"‘’s"^ = 


(252) 


(253) 


Jq 27.3 ± 1.0 days 

(Here I have simply linearly added the relative errors given for J©, about 0.0079, 


and for T©, about 0.0284, to get a conservative relative error estimate of 0.0363 for 
n©.) Multiplying hi© by the radius Rq = 6.96 x 10® m of the Sun gives Vr = 
1860 ± 70 m/s = (6.19 ± 0.22) x 10“® c. Inserting this linear surface velocity Vr and 
the escape velocity Ve above into Eq. (p51|) gives the order-of-magnitude estimate 


of Eq. (|251|) tkq as roughly 500 000 Rq, which is approximately 3.5 x 10^^ m or 

gives, which is 


2300 AU. This is about a factor of eight larger than what Eq. 
not too surprising, because of the neglect of all numerical factors and details of the 
structure of the Sun in Eq. (p5lD , and because of the high powers of the velocities 
that enter into that estimate. 

For a self-gravitating rotating fluid object (e.g., a star) which has its linear 
rotational velocity Vr/c at its surface (in units of the speed of light) greater than 
roughly the cube of the escape velocity Ve/c from its surface (again in units of the 
speed of light), as indeed is the case for the Sun, we can see from the estimate of Eq. 
(|251|) that the radius r^, where Kepler’s third law is exact for corotating circular 
orbits in the equatorial plane, is greater than the radius R of the object. If the object 
is maximally rotating, so that its rotational velocity Vr is comparable to its escape 
velocity Ve, then will be outside the object (leading to T = R^Q'^j^/{GM) > 1 
and hence a shorter orbital period at the surface than what Kepler’s third law would 
give), unless possibly the rotational and escape velocities are close to the velocities of 
light (i.e., unless possibly the object is highly relativistic, as we found was necessary 
for the polytropic models with T < 1 in the previous Section). 
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In particular, one might ask what different effective angular velocity the Sun 
would need in order that its new would then coincide with the solar radius Rq. 
To calculate this, instead of using Eq. (|251|) , use the precise formula (|25CI|) (precise 
to the extent that it gives, as it does for the Sun, an tkq that is in the far-field 
region where the field is both weak and is entirely dominated by the monopole and 
quadrupole contributions). Suppose that the effective moment of inertia of the Sun 
about its axis. 


le = 


J, 


2To 

^2 


^ = (7.12 ± 0.31) X 10^® kg m^ = (cVG')(5.29 ± 0.23) x 10^^ m® 


= (cVG) (1.25 ± 0.06) X 10^26 £3^ ^ 


3.4 X 10^ g©. 


(254) 


would stay constant as its angular velocity were changed. More precisely, assume 
that the angular momentum of the Sun would be linearly proportional to its effective 
angular velocity, and that the quadrupole moment of the Sun would be proportional 
to the square of this angular velocity, with the proportionality constants staying 
fixed. Then in order to get tk = 7?©, one would need to change the Sun’s angular 
velocity from fl© to 


n = 


Rcr 






VLq — 


3c^Q 


= (1.088 ±0.061) X 10 


-®s-' = 


27r 


© 


18.3 ± 1.0 yr 
(255) 

In other words, if the Sun were rotating with a period of greater than about 18.3 
years, then the linear term in the angular velocity (the relativistic effect linear in 
the angular momentum J) would dominate over the quadratic term in the angular 
velocity (the Newtonian quadrupole effect) at all radii outside the Sun, and so the 
corotating circular orbital period would be slightly increased everywhere outside 
such a slowly rotating Sun. 

Of course, the numerical result of Eq. (|250|) for the actual radius tkq, where 
Kepler’s third law would be exact for a corotating test body in a circular equatorial 
orbit in the gravitational held of the Sun as it is actually rotating, is entirely hypo¬ 
thetical, since the planets would exert perturbations on the orbital period far larger 
than those of APj and APj^ of Eq. (P46|) . For orbits at ~ 280 AU or greater, 
one might suppose that a reasonable estimate for some sort of averaged corotating 
period in the equatorial plane would be to use Eq. (p46|) but with the solar mass 


M©, angular momentum J©, and quadrupole moment g© replaced by the analogous 
quantities Mss, Jss, and Qss for the entire Solar System. 

Combining the data in with that in and in [Q gives directly 


GMss/c^ = 1.001 346 GMq/c^ = 1478.612 m = 0.9138 x 10®® i_ 


and 


Jss = 3.148 X 10^® kg m^ s“® = 2.985 x lO'*' h = (165.7 ± 1.3) J©, 


77 


(256) 


(257) 
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which when converted to length and area units gives 


GJss/c^ = 7.796 X 10^ = (8 829 ni)^ = 30.10 mi^ 

= 7 800 hectares = 19 260 acres. (258) 


From these data one can readily calculate that the Kerr rotational length pa¬ 
rameter a, and the corresponding dimensionless rotation parameter a, take on the 
values for the entire Solar System of 


ass = 


Jss 

Mssc 


= 52.72 km = 3.263 x lO^"" ip = (165.5 ± 1.3) Oq, 


«ss = 


assc 


cJss 

GUIs 


= 35.66 = (165.3 ± 1.3) a®. 


(259) 


(260) 


GMss 

The fact that ass > unlike the case for the Sun, means that the Solar System 
would have to give up angular momentum (in fact, give up more than 97% of its 
angular momentum) before it could possibly become a black hole. 

By adding up the time-averaged quadrupole moments of each planet and the 
Sun, around their common center of mass, from the data on the planetary masses 
(including their moons) and the semimajor axes and eccentricities of their orbits in 
], I obtained a quadrupole moment for the Solar System of 


Qss = 2.576 X 10®^ kg m^ = 1.23 x 10^° Q®, (261) 


of which about 40% came from Neptune, 23% came from Saturn, 22% came from 
Jupiter, 14% came from Uranus, 0.67% came from Pluto, 0.0026% came from Earth, 
0.0011% came from Venus, 0.00066% came from Mars, and 0.000023% came from 
Mercury. There is a positive error in my estimate from neglecting the fact that 
the orbits are not all in the same plane, and a negative error from neglecting the 
quadrupole moment contributions of the asteroids and comets, but I have not at¬ 
tempted to estimate these errors. Almost certainly not all of the four digits given 
above are correct, but I have given them just to show the answer I got for the planets 
and Sun if their orbits were coplanar. 

Converting the quadrupole moment of the Solar System to length units gives 

GQss/c^ = 1-91 X 10^^ m^ = (124 000 km)^ = 4.53 x 10^^® 4, (262) 


The effective quadrupole radius of the Solar System is then 

tqss = \J‘2Qss/Mss = 111 000 tqq = 5.09 x 10^° m = 0.340 AU = 3.15 x 10^® ip, 

(263) 

meaning that one would get the same quadrupole moment if one placed all the mass 
(of which 99.8656% comes from the Sun) into a ring at radius = 0.34 AU, 
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about 88% of the semimajor axis of the orbit of Mercury [^. Incidentally, for such 
a ring to give the angular momentum of the Solar System, it would have to rotate 
around with a period of 


Pss = = 1.03 X 10® s = 32.6 yr. 

Jss 


(264) 


The corresponding period for the 460 km solar-mass ring that gives the quadrupole 
moment of the Sun in Eq. ( 235 ) is about 14 seconds to give the solar angular 
momentum given in Eq. 


As we did for the Sun in Eqs. (|240|) , we can calculate that for the Solar System 

Qss 


Jls/Mssc^ 


= 4.65 X 10 


11 


(265) 


so the quadrupole moment of the Solar System is about 465 billion times larger than 
that of a Kerr metric with the same mass and angular momentum. 

Now if we insert these data for the Solar System in place of the corresponding 
data for the Sun alone in Eqs. (|246|) - ( P49|) , we get some sort of averaged deviations 
from Kepler’s third law for very distant orbits around the Solar System as follows: 


P = 


2tt 


n 


+K 


= Pkss + ^Pjss + ^Pqss + O(a-P^) 
271338 


= 271\ 




+ 


GMss Mssc^ 2^GMlsR 


+ 0(i?-3/2), ^266) 


^Kss 271 \ 




= 31536 986 f—) 

GMss VAU/ 


3/2 


S, 


(267) 


AF 


27rJ5s 


= 27Tass/c = (331 km)/c 


MssrP 

= 1.105 X 10"^ s = 1.105 ms 


(268) 


due to the linear effect of the angular momentum of the Solar System, and 
^t^Qss 


^Pqss - 


' AII\ 

= -1.367x 10“ (-) 

2jGMlsR ^ ^ 2 


«A2VN._15,82 (269) 


due to the quadrupole moment of the Solar System, e.g., about -2.5 days for the 
orbit of Pluto at 39.481686 77 AU [|^, assuming that the effect of the planets is 
merely to provide mass, angular momentum, and quadrupole moment for the Solar 
System. This is actually not a very good approximation for Pluto, since the formula 
above gives an orbital period of about 90 550 days for Pluto, whereas the sidereal 
period is actually 90465 days |^, about 85 days shorter. 
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One can then calculate that the effect of the quadrupole moment dominates over 
that linear in the angular momentum for the Solar System for 

r < tkss = ^ = 2-29 X lO^^m = 1.53 x 10^^ AU = 7.42 x 10®Mpc. (270) 

Since tkq dehned by Eq. ( 115U|) for the Sun alone is larger than the Solar System, 
and since the radius r^ss which Kepler’s third law would be exact for a test 
body orbiting the entire Solar System in otherwise flat spacetime is far larger than 
the presently observable universe, we can conclude that for realistic orbits (i.e., at 
orbital distances less than that to the next nearest star), the quadrupole moment 
of the Sun or of the planets always dominates over the linear effect of the angular 
momentum in changing the period of circular orbits from the value given by Kepler’s 
third law, so that the period is always smaller than that given by Kepler’s third law. 
(This is under the approximation that the effect of the planets is merely to give a 
quadrupole moment to the Solar System and ignores more complicated effects when 
the orbital radii, and hence periods, of the planets are non-negligible fractions of 
the radius and period of the orbit of the test body.) 

However, we saw in the previous Section that for certain relativistic polytropic 
star models with polytropic index not too large (not too soft an equation of state), 
it is possible to have the linear effect of the angular momentum dominate over that 
of the quadrupole moment, even at the surface of the star, so that the corotating 
Keplerian orbits at the equatorial surface of the star can have a longer period (as 
seen by a nonrotating observer at inhnity) than what Kepler’s third law would give. 
As discussed above, this effect is related to the fact that the angular momentum of a 
source generally causes the Stationary Congruence Accelerating Maximally (SCAM) 
to be counterrotating, so that for a given positive orbital angular velocity relative 
to the SCAM, the angular velocity relative to infinity is less, giving rise to a longer 
period. 
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